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Abstract 

Based on a general variational principle, Noether's theorem is re- 
visited. It is shown that the so called pseudotensor problem of the 
gravitational energy-momentum is a result of mis-reading Noether's 
theorem, and in fact, all the Noether's conserved quantities are scalars. 
As a side product, a generalized Hamilton- Jacobi equation for the 
Hamilton's principal functional is obtained. 

1 Introduction 

Ever since Einstein introduced the gravitational energy-momentum complex 
to keep the law of conservation of energy-momentum alive in general relativ- 
ity (GR) in 1918, the so-called pseudotensor problem of various gravitational 
energy- momentum complexes has been troubling relativitists [1]. It is widely 
accepted that the proof of the positivity of the total gravitational energy both 
in spacial and null infinity is one of the greatest achievements in classical GR 
in the last 30 years. This success inspired efforts to find quasi-local gravita- 
tional energy-momentum. Non-locality is definitely not acceptable; however, 
finding quasi-local gravitational quantities has proven to be surprisingly dif- 
ficult [2]. Experimentally, the gravitational wave carrying energy-momentum 
has not been tested [3]. The aim of this article is to explore the pseudotensor 
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problem of the gravitational energy-momentum, in particular, the difficulty 
of non-locality of the gravitational energy-momentum. The whole argument 
is based on the variational principle of dynamics. It is pointed out that 
the so-called pseudotensor problem of the gravitational energy-momentum 
is only a result of mis-reading Noether's theorem, and mistaking different 
geometric physical objects as one and the same; hence the non-locality dif- 
ficulty does not exist at all. In fact, all Noether's conserved quantities are 
scalars. As a side product, a generalized Hamilton- Jacobi equation for the 
Hamilton's principal functional is obtained. 

The history of modern physics has proven that the variational principle 
approach to dynamics is not only an alternative and equivalent version to the 
naive, intuitive approach, but also yields deeper insights into the underlying 
physics. For instance, it is hard to imagine that the statistical mechanics 
could have been established without using the concept of phase space, and 
that the quantum mechanics could have been established without using the 
concept of Hamiltonian. Therefore, we will base our argument on a gen- 
eral variational principle for the classical field. It might be for the similar 
consideration that soon after Einstein proposed his general theory of rela- 
tivity, Hilbert made the first succesful attempt to get Einstein's equation 
by using the least action principle. The Lagrangian being used for vacuum 
Einstein's equation, (16^(7) is the only independent scalar constructed 
in terms of the metric field and its derivatives of no higher than the second 
order. However, because the Ricci scalar curvature R contains the second 
order derivatives of the metric field g /Jil/ (x), which is now the dynamic vari- 
able along with the matter field, the least action priciple for Lagrangians 
containing only the field quantity and its first order derivatives does not lead 
to Einstein's field equation. The generally accepted solution to this diffi- 
culty is adding the Gibbons- Hawking boundary term to the Hilbert action 
and keeping the least action principle unchanged [4]. But there is another 
solution to this difficulty, which was initiated by Hilbert and will be adopted 
in the present paper. The least action principle will be restated and the 
Hilbert action will still be used for the vacuum Einstein's equation. In order 
to show this is proper and natural, we will consider the variational principle 
for classical fields in 2)-dimensional spacetime with a Lagrangian con- 
taining the field quantity and its derivatives of up to the l)-th order. 
In our opinion, non-local interaction is not acceptable, so we assume that the 
Lagrangian does not contain the integration of the field quantity. 

There have been varied versions of variational principle and Noether's 
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theorem in the literature [5], and different notations have been used by dif- 
ferent authors. For the readers' convenience, and for the consistency of the 
reasoning, we start with presenting a general variational principle for classi- 
cal fields in 2)-dimensional spacetime with a Lagrangian containing the 
field quantity and its derivatives of up to the iV(> l)-th order. In section 2, 
the Lagrangian formalism is presented, and a generalized Hamilton- Jacobi 
equation for Hamilton's principal functional is obtained. In section 3, the 
Noether's theorem is rederived. Then, in section 4, these general results are 
applied to the specific case of general relativity, especially Noether's theorem 
is applied to get quite a few conservation laws in GR. It is noted that most 
of the conserved quantities obtained here are not tensors like various gravi- 
tational energy-momentum complexes. Then Noether's theorem is revisited 
in section 5, and it is shown that the so-called pseudotensor problem of var- 
ied gravitational energy-momentum complexes is a result of taking different 
geometrical physical objects as one and the same; as a matter of fact, all 
the Noether's conserved quantities are scalars; and the non-locality difficulty 
does not really exist. 

These results will be used to explore the energy-momentum conservation 
and the gravitational energy-momentum in GR in a later paper. 

2 A general variational principle for classical 
fields 

2.1 Lagrangian formulation 

First, we present a useful mathematical formula for the variational principle, 
which does not rely on physics. 

Suppose {$ B : W 1 — > R (or C) | B = 1, 2, . . . , /} (n > 2) are smooth 
functions, and function L = <9$(:r), <9 2 $(:r), . . . ,d N &(x)) is smooth 

with respect to all its arguments. It is easy to show just by using Leibniz's 
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rule that (the Einstein convention is used for coordinate indices)Q 

f N nr 



B .., v . „ dd Xl ---d Xx &B(x) 

f n B f N-l 

E D-'^ • • ■ 9 ^w-^m + * E E 

B=l X=0 Ai Ax ^ v 7 B=l Z=0 

N-l-Z or 



••*,•„(*> E '-"' 8 '.-^ w ,,.. 8 , r4 ,... 4l ^M 
^^)D-i)\A ^,„ M , ( „ 

/ N-l 

+ dxJ2J2 5d ^-" d» z ®B{x)K BXv ^{x) (1) 



B=l Z=0 



when $b(^) 1 — > $b(x) = $b(x) + 5§ B (x),VB = 1, ...,/. Consider the 
functional F of the following form 

= [ d n xL(x,<$>{x),d<$>{x),d?<$>{x),...,d N ${x)), (2) 

where f2 is an open subset with a compact closure of lR n . When the arguments 
$b(^) V-B = 1, . . . , / change slightly, the variation of functional F is 

f N dL 



5F[$]= I d n xJ2J2 6d ^--- d ^ Bi 

J n B=1X=0 



X 



ddx 1 ---dx x $ B (x) 



1 In order to avoid the indefiniteness of derivatives such as gg^y L, and in order to 
keep the formulaes neat, it is assumed in the present paper without loss of generality that 

L{x, $(x),d$(x), d 2 1>(x), d N $(x)) 
= : L(x a ,$ B (x),d ai $ B (x),d ai d a2 $ B (x), . . . ,d ai d a2 ■ ■ ■ d aN $ B (x)) 

= L(x a , $ B (x),d a <$> B (x), — ^ PeSa d apw d ap{2) <I> B (x), 



We will pretend that all the cross derivatives are independent variables of L, when calcu- 

d 

dd 1 d 2 'S>E 



lating aa a & L, etc. See Appendix A for the details. 
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, f N-l 

+ / ^xdxY.Y, 69 ^ - ■ ■ d ^^B{x)K B ^- Uz {x) 

J Q B=l Z=0 



1 V 



/ AT-l 

■"*(*) (3) 



, / N-l 

+ / ds A (x) ^ Sd »i ■ ■ ■ dv z <5> B (x)K BX ^ 

J 80, B=1 Z=Q 



This can be easily obtained by using eqn.(l) and the Stokes theorem. The 
derivatives of functional (2) is defined as follows. 

Definition 1 If the change of the functional (2) can be expressed as 

r f 

F[$ + 6$] - F[$] = d n x^2S^ B (x)D B [^,x]+o[S^} (4) 

J Q B=l 

where D B [&,x] is a functional of $ varying with x, and o[5Q] is a higher 
order infinitesimal of 5$, when | B — 1, 2, . . . , /} change slightly while 
the boundary values of $, <9$, . . . , d N ~ l Q are kept fixed, then F is called dif- 
ferentialble at $ 7 and D B [§,x] is called the derivative of functional F with 
respect to <&b at <3> and point x, and denoted by 

Let us now apply the general formula (3) to the action functional of 
classical field {$# : M — > R (or C) | B = 1, 2, . . . , /} 

A[$]= f d n xL(x,$(x),d<S>(x),d 2 <S>(x),...,d N $(x)) (6) 

Jx(Q) 

where Q is an open subset with a compact closure of the spacetime manifold 
M, x(Q) C M. n is the image of Q C M under the coordinate mapping x : 
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M— > M. n and L is the Lagrangian of the field. We get the difference between 
the action functionals over Q of two possible movements (two paths allowed 
by the constraints) close to each other 



r f N ar 

6A[*]= / ers£*Ms)£(-l)^---^ ^T^y 



f N-l 

■■•"(*) 



+ / d n xd x J2J2 Sd ^--- dv z <5> B (x)K BX ^ 

Jx(P) B=1 Z=Q 

r f N 8T 

= d n x ^ B { X )Y J {-i) x d^---d Xx - — —— 

Jx(Q) tZ Z~ a O0 Xl - ■ ■0 Xx <Pb{X) 



B=l X=0 
r f N-l 

+ / ds x ( X )j2J2 6d ^--- d ^ B ^ KBXvi " Mz ^ w 

Jx(dQ) B=1 Z=Q 

Equation (7) suggests that for all AT ^ 1, n ^ 2 the least action 
principle read as follows. 

For any spacetime region Q, among all possible movements in Q with the 
same boundary condition 

8<Z>\ dn = 0, 5d<$>\ dn = 0, . . . , 6d N ^^\ dn = 0, (8) 

the real movement (the path allowed by physical laws) corresponds to the 
stationary value of the action over Q. 

Combining eqns.(7), (8), one obtains the field equation (Euler-Lagrange 
equation) satisfied by the real movements 

This Lagrangian formalism applies to Newtonian dynamics, dynamics in 
special relativity (SR) and dynamics in general relativity. Specific covariance 
is the heritage from the Lagrangian being used. It is not the heritage from 
this general formalism. 



6 



2.2 Hamilton's principal functional and Hamilton-Jacobi's 
equation 

Let us consider the difference between actions over spacetime region Q of two 
real movements close to each other. Using eqns.(7) and (9), one gets, for real 
movements 

, / N-l 

6A[$] = / ds x (x) E 5d »i • ■■dv z $ B (x)K BX ^(x) (10) 

From eqn.(10), one sees that the action over a spacetime region Q of a 
real movement is determined by the closed hyper-surface dQ, and 
d^\dn,..., d N ~ 1 &\df l . It will be called the generalized Hamilton's princi- 
pal functional and denoted by 

S = S[dn, Siangan, ... , d^ 1 ^} (11) 
Re- write eqn.(10) as 

n f N-l 

ss= ds x (x)Y,J2 6d ^--- d ^ B ^ KBX,/1 ^ z ^ ( 12 ) 

Jx(dQ) B=1 Z=Q 

Note that when ^slon is given, only one of the n derivatives d\^s\dn 
(A = 0, 1, . . . , n — 1) is independent; when d\§B\an is given, only one of 
the n derivatives d^dx^slan (// = 0, 1, . . . , n — 1) is independent; and so on. 
Thus for a given suffix B, only iV items ( Qsldn, one of dx^slan, ■ ■ ■ , one of 
<9aj ■ • • d\ N _ 1 ^ B \dn (Aj = 0, 1, . . . , n - 1) ) are independent. 

In order to formulate the generalized Hamilton-Jacobi's equation, one 
needs a new type of functional derivative [6]. 

Definition 2 Let I] be a closed hypersurface in W l , \I> a function defined on 
R n , and F = F[E, a functional o/E and The functional derivatives 
are defined as follows. If the variation of F can be written as 

<JF[E,*| E ] = y ^ A (x){F[E,^| E ,x]^(x) + Z[E,^| s ,x] A ^(x)} (13) 

where F[E, x] A and Z[E, x] A are junctionals of E and which 
vary with x, when E i — >■ E ; and \1/ i — >■ \l/ 7 £aen taey are called the functional 
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derivative of F with respect to E M (x) and ^{x), and denoted by 

respectively. 

Hence we have 



5F[E,*\v]= J ds x (x) 



SW(x) + 



5^(x)J ' y 1 \S^(x) 
The hypersurface E is given by the parameter equation 



5V(x) 



x r 



The <5E M (:r) and S^(x) in eqn.(13) are respectively 

Now, from eqns.(12) and (15) we get [dQ in (12) is E in (15)] 



(15) 



(16) 



(17) 



5S 



5d Vl ■ ■■d Uz $ B {x) 



K BXvi ~ Vz (x), 
0,...,N-1,B = 1,...J. 



(18) 



Following the evolution of one real movement and observing the change of 
its action with E, 



5S= I ds x (x)L5 x JJ:' T (x) 
'x(dn) 

f N-l 

+ EE 



STTix) 



B=0 Z=0 



SS 



5d vi ■ ■■d Uz ^ B (x) 



d a d Ul ---d Uz <S> B (x)5Z°{x)} (19) 
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We get the generalized Hamilton- Jacobi's equation. 

A / JV-1 

+ l LS ° "EE d ^ ■ ■ ' d^B{x)K BX ^{x)\ = (20) 

B=0 Z=0 

It will be shown in the next section that the expression in the bracket is 
just the canonical energy- momentum t*(x), hence we have 

For a dynamic system with finite degrees of freedom, we have the Hamilton's 
principal function S^q 1 , . . . ,q h ,t), and the Hamilton-Jacobi's equation is a 
partial differential equation || + H(q,^-,t) = 0. While for the classical 
field, we have the Hamilton's principal functional, and the Hamilton-Jacobi's 
equation is a set of differential equations for functional. It plays an important 
role in canonical quantization for fields. 




3 Noether's theorem 
3.1 Re-deriving the theorem 

There have been varied versions of Noether's theorem, and different notations 
have been used. To avoid confusion, here we redrive Noether's theorem for 
classical fields in n-dimensional spacetime with a Lagrangian containing field 
quantity and its derivatives of up to the iV-th order, in terms of coordinate 
language. 

Theorem 3 If the action of classical fields over every spacetime region Q 
with a compact closure 

A[$]= [ d n xL(x,$(x),d<5>(x),d 2 $(x),...,d N <5>(x)) 
remains unchanged under the following r— parameter family of infinitesimal 
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transformations of coordinates and field quantities 

x x i — > x x = x x + Sx x , 
Sx x = : S x x (x, e 1 ,..., e r ), |e*| < 1, 5x A (x, 0, . . . , 0) = 
$ B (x) i — >■ = $ B (x) + 8$ B (x), 

5<Z> B (x) =: 5$b(x, e\ . . . , e r ), ^ < 1, 0, . . . , 0) = (22) 

£/ien there exist r conservation laws. 

Proof. The small change of field quantity 6<&b{x) can be devided into two 
parts, the part due to the small change of its function form and the part due 
to the small change of its arguments respectively. 

£$ B ( x ) = M B (x) + 5x a d a <Z> B (x) (23) 

Similarily the small change of derivatives of field quantity S[d Xl d X2 ■ ■ ■ 8 Xx ^b{x) 
can be written as 

5[d Xl d X2 ■ ■ ■ d Xx <5> B (x)] = d Xl d X2 ■ ■ ■ d Xx 5$ B {x) 

+6x°d a d Xl d X2 ---d Xx <S> B (x) (24) 

The variation of the action can be devided into two parts. One is due to the 
small change of the integration domain x(Q) i — > x(Q) in IR ra , and the other 
is due to the small change of the integrand 

L(x, $(x), d$(x), d 2 <S>(x), d N $(x)) 
i— >■ L(x, §(x),d$(x),d 2 $(x), . . . , d N $(x)) 



5A{$] = [ ds x (x)5x x L 
Jx(dn) 



/ N 



+ / dn * E E dd d~$TJx) d ^ ' ' ' 

/ N 



dL 



[ ^^^(x)^(-l)^-.-^— 



B=l X=0 

f N-l 

I , r - \ - 1 



/ d n xd x [5x (T 5 x L + J2J2 d ^--- d -z^B(x)K BX ^- Vz (x)} (25) 

Jx(P) B =l Z=0 
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The first integrand at rhs vanishes for real movement, hence the second 
integral does too. One gets the following continuity equation due to the 
arbitrariness of Q. 

f N-l 

d x [5x°5 x L + £5> 1 -"^ (S*b(x) - 8x°d a <S> B {x)) K BXvi ~ Mz (x)\ = 

B=l Z=0 

or 

, f N-l 

/ ds x (x){5x^5 x a L + J2Yl 

[d vi • • • 3„ z {5<S> B (x) - Safd ff ^ B (x)) K BX ^(x)}} = (26) 

Noting that both 5x a and 5<& B (x) depend on r real parameters, one gets from 
eqn.(26) r conservation laws. 

B=l Z=0 ^ ' 

dSx a 

+ —\ e=0 5 x L} = 0, Va = l,...,r (27) 



3.2 Conservation law due to "coordinate shift" invari- 
ance 

In this subsection, we restrict our discussion to those cases, such that the 
Lagrangian does not manifestly contain coordinates and the action of classical 
field over any spacetime region 

A[$]= f d n xL^(x),d^(x),d 2 ^(x),...,d N ^(x)) (28) 

Jx(Sl) 

remains unchanged under the following "coordinate shift". 

5x a = e a , 5$ B (x) = 0. (29) 

And eqn.(27) reads 

d x r x (x) = 0, (30) 
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where 



/ N-l 

r x a (x) = ^"££ d ^ ■ ■ ■ 9 ^b{x) K B ^^(x) (31) 

B=l Z=0 

is usually called canonical energy-momentum. 

Comparing eqn.(20) and (31), we get eqn.(21), the generalized Hamilton- 
Jacobi equation. 

The formalism presented so far is good for any classical field in n-dimensional 
spacetime with a Lagrangian containing field quantities and their derivatives 
of up to the iV-th order, no matter it is Galileo covariant, Lorentz covari- 
ant, general covariant or without any covariance. Therefore, we can not say 
eqn.(21) is covariant even though it looks like so. The covariance of the for- 
malism is the heritage from the Lagrangian being used, not the heritage from 
the formalism itself. 



4 Variational principle approach to general 
relativity 

Let us apply the general results obtained above to the classic fields in GR. 
We will consider the case where the matter field is a (1, l)-tensor field. The 
results can be readily generalized to any (r, s)— tensor matter field. For the 
dynamic system, (1, l)-tensor matter field u°(x) plus the metric field g Q( g(x), 
the Lagrangian and the action over any spacetime region Q are 



L(g(x), dg(x), d 2 g(x), u(x), du(x)) 

^-\g{x)\[C{g{x),u{x),Vu{x)) + R\ = L M + L G , (32) 



16itG 
Jx(n) 

+ j^Q R ] = a m[9, y] + Mg\ (33) 

where R is the Ricci scalar curvature, C(g(x), u(x), Vm(x)) is the Lagrangian 
for matter field obtained from the Lagrangian in special relativity 
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C{j], u(x), du(x)) by replacing the Minkowski metric 77 = diag(—l, 1,1,1) 
with g a p(x), and replacing the partial derivatives d^u®(x) with the covariant 
derivatives V m m|(x). 

4.1 Einstein's field equation 

The Euler-Lagrange equation, Eqn.(9) now reads 

5A\q,u] i — . . . dC „ dC . , , 



(35) 



where T a @ is the energy-momentum tensor of matter field, which is a sym- 
metrical (2,0)-tensor. 



0V A uf(aO fy^x) 1 dV ,u\{x) dd x g aP {x) 

_ 2d dC dV^x) _ 2 d£ g gV^jgx) 

A 5V^(a;) dd x g a p{x) <9V M u|(x) A d<9 A # a/3 (:r) 

= T^( M (x),a M (x),^(x),^(x)) (36) 

All the pre-GR dynamics discuss how matter behaves in spacetime with a 
given metric, that is, Newton's absolute space-time or Minkowski space. All 
of them can not explain why there are inertial reference systems and non- 
inertial reference systems, and why the inertial mass equals the gravitational 
mass. G.R. is unique. It discusses how the matter motion determines the 
bending of the spacetime and how the matter behaves in a curved spacetime. 
Hence the dynamical variables are both the matter field u^{x) and the met- 
ric field of spacetime g a p(x). It means that they are both to be determined 
simultaneously by solving the equation of motion. It does not necessarily 
imply that the metric field of spacetime g a p{x) is matter (in its philosophic 
context, in the most general meaning of the word) like, say, the electromag- 
netic field, the spinor filed, etc. It does not necessarily imply that the metric 
field of spacetime g a p(x) carries energy- momentum. 
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4.2 Noether's theorem for classical field in GR 

The Noether's conservation law, or the continuity equation (26) , now reads 

{V-\9i^)\J K [u(x),du(x),g(x),dg(x),d 2 g(x)]5x,8u(x),8g(x),d5g(x)]} =0 

or 

ds\(x){y/-\g(x)\ 



lx(dQ) 

J K [u{x),du(x),g{x),dg{x),d 2 g{x);5x,5u{x),5g{x),d5g{x)]} = (37) 
where 

J K \u(x),du(x),g(x),dg(x),d 2 g( y x);5x,5u(x),5g(x),d5g(x)~\ =: J£ 

1 m^+i^J^-o, m 



16ttG p " dd K g a p(x) " dd K d^g a p(x) 
OR - 5 dR 



- T U X ) ^ p „ .,J 5 9*e(*) + fl : _ ( ^ S gaP {x)} (38) 



It is worth noting that the form of function J K is independent of coordinate 
systems, and the arguments of function J K are u(x), du(x),g(x), dg(x), d 2 g(x); 
and Sx,5u(x),5g(x),d5g(x). 

4.3 Conservation laws in GR 

4.3.1 Conservation law due to "coordinate shift" invariance 

Action (33) remains unchanged under the following "coordinate shifts". 

5x p = e p , 8ul(x) = 0, 8g aP (x) = 0. (39) 
In this case, eqn.(26) reads 
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dx[V Z \dWK(x)} = 0, (40) 

where 



T p( x ) = T K p (u(x),du(x),g(x),dg(x),d 2 g(x)) 

= n 8C a , , d p u d Jx) - C5 K n + - L dC , , g 6a (x)vPAx) 
dW K u\{x) p ^ J p 2 [ dW K u\{xY y ! ^ ! 

+ RX7 9C e( , 9 ea (x)u«(x) - ^ r <f"(x)ul(x) 

- n^(x)uUx) - —^—g^( x )u e Ax) 

dC 1 

U Q Til' / \ \ Q / \ 

" ^ddjtfapix) ~ ^ {X) dd K d,g aP {x) )dpga ^ X) 

+ a^M BMx) - RS;] (41) 

is usually called canonical energy- momentum tensor. It is worth noting that 
the arguments of function T^ ) ( y u(x),du( y x),g(x),dg(x),d 2 g(x)) are all field 
quantities and their derivatives. Notice that "coordinate shift" eqn.(39) 
is not an invariant concept under general coordinate transformation. This 
explains why r A (x) is not a tensor under general coordinate transformation. 
We will get back to this problem later. 

4.3.2 Conservation law due to "4-dimensional rotation" invari- 
ance 

The action (33) remains unchanged under infinitesimal "4-dimensional ro- 
tations" (Lorentz transformtions), which form a 6-parameter family of in- 
finitesimal symmetry transformations 

^ ^ ^ = L(X) L , = ^ + A m \ K \ « 1; ^ A A = V) A ^ = 
5x A = A A ^ = i(r7 A "x CT - ?7 A,T x p )A p(7 , (42) 
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6u°(x) = A^(x) - Ay v (x) 

'rfi p u l . 

2' 



K p J-[rf»ul{x) - V e °u P (x) - S^ p u e v (x) + 5 p rj^u e v (x)} : 



5g a p{x) = A^l-S^g^x) 

+ Kv^gA*) ~ W P 9cu{x) + 5^g au (x)] 
In this case eqn.(26) reads 

d 

dx K 

where 



{^\g\x)\M^} = Q 



2M Kpa (u(x),du(x),g(x), dg(x), d 2 g{x)) 

f)C 

= CS^x" - rfxr) + — A W P ul{x) - rf°if(x) 

Ov K U^ [X ) 

- S^u^x) + 8lrf°u 6 v {x)\ - {r] Xp x c - r] Xc x p )d x u\{x)} 
+ - L fl/ J a {x)uHx) + - ftl - g ea (x)u K Ax) 
()C ^k ( ^B(^ 9C j-« (a .\,fi, 



{[- (*) + ^rTg^x) - 5X P 9au(x) + ^^(x)] 

- - ^x p )d x g a Ax)} + ^ " V Xa x p ) 

+ (_?« a _V (x) ^ ) 

dd K g af3 (x) * dd K d^g a/3 (x) ">* dd K d^g a/3 (x) 

(l-Kv» p g^(x) + Sfrrgtf(x) - SX P 9au(x) + S p ^g av (x)} 

BR 

- - ^)^(x)) + 

+ SirTg^x) - 5^ p g au {x) + 8 p p rf a g av (x)} 
-(r^x* - v x,T x p )d x g a p(x))} 
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4.3.3 Conservation law due to "4-dimensional pure deformation" 
invariance 

The action (33) remains unchanged under infinitesimal "4-dimensional pure 
deformations", which form a 6-parameter family of infinitesimal symmetry 
transformations 

v^t = V> V = A ^ ( 47 ) 

Sx x = A x x p = \{r] Xp x a + r] Xa x p )A p(T , (48) 
6u°(x) = Aj«?(x) - AJ«5(x) 
= + " ^V P <W - 5 P r)^u e v (x)]A pa , (49) 

+ S p a rr9rf(x) + V p g ai/ {x) + ^r/^^,(x)]A, (7 (50) 
In this case eqn.(26) reads 

d 



dx K 



{^\gJxJ\N^} = (51) 



where 



2N Kpa (u(x), du(x), g(x), dg(x), d 2 g(x)) 

fir 

C5 K x ( V Xp x° + ^x p ) + A W P ul{x) + rfu\{x) 

O V K U^ yX ) 

-5lrf p u e v {x) - 5 p r]^u v (x)} - ( V xp x a + rf* x p )d x u\{x)} 
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OR dR * dR 



dd K g a/3 (x) dd K d^g aP (x) Vfi dd K d^g a ^x) 
+ 5 p ^g av {x) + (77^ + 

+ 8 P pn va goa ,{x) + ( V Xp x° + v Xa x p )d x g a p(x)}} (52) 

4.3.4 Conservation law due to "scaling" invariance 

The action (33) remains unchanged under infinitesimal scaling transforma- 
tions, which form a 1-parameter family of infinitesimal symmetry transfor- 
mations 

x x \ — >x x = e e x x , |e| < 1, fe A = ex A , VA = 0,1,2,3 (53) 
= 0, 8g a p(x) = -2eg a/3 (x) (54) 

In this case eqn.(26) reads 

d 



dx K 



{y/-Wx)\S*} = (55) 



S K (u(x), du(x), g(x), dg(x), d 2 g(x)) 



r)/" ar 



+ _}—i R s« x p _ r ^ a ^ r fx) 9jR 



IGttG 1 ' [ ^^ Q/3 (x) » ddjtfapix) upK >dd K d,g aP {xy 

dR 

[2g a p(x) + x p d p g a p{x)\ - — - - dpftgapix) + x p d p g a p(x)]} (56) 
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4.3.5 Conservation law due to "skew-scaling" invariance 

The action (33) remains unchanged under infinitesimal "skew-scaling" trans- 
formations, which form a 1-parameter family of infinitesimal symmetry trans- 
formations 

Sx° = - e V, Sx 1 = e V, 5x 2 = 0, Sx 3 = (57) 

dx 
dx 

5u° {x) = 0, 5u1(x) = -2e 1 M°(x), 5u° 2 (x) = -e^x), 5u° 3 {x) = -e^x) 
Sul(x) = 2e 1 ul(x), 5u\(x) = 0, 5u\(x) = t l u\{x), 5u\(x) = e 1 ^^) 
Su J (x) = e^^x), 5u{(x) = —e^^uKx), 5u 3 k (x) = 0, Vj, k = 2,3 
$9ao(x) = Z^gaoix), $9oi(x) = 0, Sg 02 (x) = e 1 g 02 (x), 5g 03 (x) = e 1 g 03 (x) 
$9io(x) = 0, 6g n (x) = -Z^guix), 8g 12 (x) = -e 1 g 12 (x), 8g 13 (x) = -e 1 g 13 {x) 
Sgjo(x) = e 1 g j0 (x) : 5g jl (x) = -€ l g jl (x), 5g jk (x) = 0, Vj, k = 2,3 (59) 
Substitute eqns.(57) and (59) into eqn.(26), we can get a conserved current 

^{V=5Sl^} = (60) 

Here we skip the expression of J£(u(x),du(x),g(x), dg(x), d 2 g(x)). Similarily 
we can get conserved currents J 2 and J 3 . 

Because the symmetry group of classical field in GR is an infinite di- 
mensional Lie group, we can get infinitely many conservation laws by using 
Noether's theorem. It is interesting to note that most of the conserved cur- 
rents are not tensors. For instance, r n p is not a (1,1) tensor, M Kpa is not a (3,0) 
tensor, N Kpa is not a (3,0) tensor, and S K is not a (1,0) tensor, etc. One of 
the main tasks of the present work is to show that these pseudotensor prob- 
lems (including the pseudotensor problem of various gravitational energy- 
momentum complexes) are just results from misreading Noether's theorem. 



1-e 1 








1 + e 1 

1 

1 



(58) 
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5 Noether's theorem revisited 



Let us start with the simplest example. The canonical energy-momentum 
Tp(u(x),du(x),g(x),dg(x),d 2 g(x)) in eqn.(41) is not a tensor 

T K p (u(x),du(x),g(x), dg(x), d 2 g(x)) ^ 

This can be understood as follows. The 4-parameter families of infinitesimal 
transformations of coordinates and field quantities eqn.(39) and 

5y» = e", 8u\{y) = 0, 5g af) (y) = 0. (39') 

do not correspond to the same 4-parameter family of infinitesimal diffeomor- 
phisms of spacetime M, when 

3 p e M, and < a, (5, a < 3, such that - y — ? L ^ 0, (62) 

dx a ox p 

Therefore, generally speaking, the conservation laws (continuity equations) 
corresponding to them are not equivalent to each other, and the conserved 
quantities corresponding to them are not the same geometrical physical ob- 
ject. Inequality (61) is a comparison between components of two different 
geometrical physical objects which have been mistaken one and the same. 

In general, infinitesimal transformations of coordinates and field quan- 
tities expressed in different coordinate systems with the same form (like 
eqns.(39) and (39')) do not correspond to the same infinitesimal diffeomor- 
phism of spacetime M. Conversely, the same r-parameter family of infinitesi- 
mal diffeomorphisms of spacetime M is described differently in different coor- 
dinate systems. For instance, if a 4-parameter family of infinitesimal diffeo- 
morphisms of spacetime M is described in coordinate system ( 
by eqn.(39), then it would be described in coordinate system (y° , y l , y 2 , y 3 ) 
not by eqn.(39'), but by 

<Wy) = -^(£)9M + ^-/£i)9 P ^)\ (39") 
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which generally is no longer "coordinate shift" . 

It is generally accepted that in GR, the same physical geometrical object 
is expressed in different coordinate systems by functions of field quantities, 
their derivatives, and coordinates, with the same function form. But 
this is not always true. For instance, the Noether's conserved current (38) 
J K \u{x), du(x), g(x), dg(x), d 2 g(x); 5x, 5u(x), Sg(x), dSg(x)~\ , and its counter- 
part corresponding to the same r-parameter family of infinitesimal diffeo- 
morphisms of spacetime M, but written in coordinate system (y°,y l ,y 2 ,y 3 ), 
J K [u(y),du(y),g(y),dg(y),d 2 g(y)]5y,5u(y),5g(y),d5g(y)], are not function 
of field quantities, their derivatives, and coordinates, with the same function 
form. Because J K manifestly contains 5x, 5u(x), 5g(x), dSg(x), besides field 
guantities, their derivatives and coordinates, while 5x, 5u(x), Sg(x), dSg(x) 
and 5y,5u(y),5g(y),d5g(y) are different functions of field quantities, and 
their derivatives, and coordinates. Yet they are the same physical geometri- 
cal object, as will be shown in the following. 

To prove this, let us show that continuity equations written in different 
coordinate systems, but corresponding to the same r-parameter family of 
infinitesimal diffeomorphisms of spacetime M , are equivalent to each other, 
even though the conserved currents written in different coordinate systems 
are different functions of matter field, metric field, their derivatives, and 
coordinates. 

Let {A(e 1 , . . . , e r ) =: A : M — )> M | |e*| < 1,V1 < i < r} be an r-parameter 
family of infinitesimal diffeomorphisms of spacetime manifold M onto itself, 
and A(0, . . . , 0) the identity mapping of M. For any chart (U, ip) of M, de- 
noting by ( 3 ) the corresponding coordinate system, we have 

x\ x 2 , x 3 ) ^T^" x x = x x + 5x x (x; e\ . . . , e r ) =: x x + 5x x , (63) 
where 7r A : (x°, x 1 , x 2 , x 3 ) x x is the projection operator, 

5ul(x;e\...,e r ) = : 5u\{x) = ^^u^(x) - ^^u 9 v (x), 
8g a p(x;e 1 ,...,e r ) = : 6g aP (x) = -^^g^{x) - ^^-g av {x) (64) 

and the continuity equation (37) which is re-written for simplicity as 

d 



dx 



-[yfW)\J2\ = (65) 
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Switching to coordinate system (y°, y 1 , y 2 , y 3 ), which corresponding to chart 
(V, ip), we have 

y x + 5y x (y,e\...,e r )=:y x + 8y, (63') 



{y\y\y 2 ,y 3 ) .— >• y> 



Sg a fi(y,e L ,...,e r ) 



dyV ™> Qyt V 
d5y» dby v 



dy c 



dyP 



d 
dy K 



V z \gW\Jy 



= o 



9aM (64') 



(65') 



For the same r-parameter family of infinitesimal diffeomorphisms of space- 
time manifold M onto itself {A (e 1 , . . . , e r ) =: A : M ->■ M | |e*| < 1, VI < % < r}, 
we have 



8y x = d ^°J 0( P) 5 ^ = 9 J^8x 



(66) 



dx K dx K 

From this relation it is easy to see that Sx x (x; e 1 , . . . , e r ) and Sy x (y; e 1 , . . . , e r ), 
generally speaking, are functions of different forms. However, we will show 
that the continuity eqns. (65) and (65') are equivalent to each other. J£ and 
Jy are respectively the components in coordinate systems (x°, x 1 , x 2 , x 3 ) and 
(y°, y 1 , y 2 , y 3 ) of the same vector field. 



Theorem 4 



dx 



(67) 



Proof. The lhs of eqn.(65) can be written as 



dy 



dx 



JK 



dy x d 
dx K dy x 

9 f V- 



dy> 



dy x 



dy 



\9 



d 2 y x dx a 
dx K dx a dy x 
d 



JK 

dx K x 



dy 



dx 
dy 



d 



dx K 



dy 



dx 



dy> 



dx 
dy x 
dx K 



dx 



d dy x 
dy x dx R 



dy 



dx 



JK 





dy 




dx 



(68) 
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here an identity on the Jacobian 



d 

dx K 



dy 



dx 



d 2 y x dx a 



dy 



dx 



d dy x 
dy x dx K 



dx K dx a dy x 

has been used. Hence the proof is reduced to proving 

dy x 



dy_ 

dx 



- J K 

dx- x 



j: 



X 



u^{x) 



(69) 



(70) 



In fact, it's easy to show that 5x K is a vector, 5u®{ 
5u®(x) — Sx p d p u^(x) is a (1, l)-tensor, and Sg a/ 3(x) = g a p(x) — g a /3( x ) = 
5g a p(x) — 5x p d p g aj 3(x) is a (0,2)-tensor. Hence the terms in the first brace 
of J£ eqn.(38) are vectors, and the first term in the second brace is a vector 
too. The rest of the terms in the second brace are not vectors individually. 
However, their sum is a vector. This can be proven straightforwardly, though 
tediously. (See Appendix B) ■ 



Conclusion 5 Now we have proven: 



(i) It is not the Noether's conservation laws written in different coordinate 
systems corresponding to the infinitesimal coordinate transformations with 
the same form that are equivalent to one another; rather, it is the Noether's 
conservation laws corresponding to the same family of infinitesimal diffeo- 
morphisms of spacetime onto itself that are equivalent to one another. 

(ii) All the Noether's conservative currents are vector fields (eqn.(70)), 
which should be the density, and current density of some scalar. Therefore, 

(Hi) All the Noether's conserved quatities are scalars. 

It is worth noting that the form of function J K eqn.(38) as a function 
of u(x),du(x), g(x),dg(x),d 2 g(x), Sx,Su(x),Sg(x),dSg(x) is independent of 
coordinate systems. However, for a given family of infinitesimal diffeomor- 
phisms of spacetime, the function forms of 5x(x, e 1 , . . . , e r ), 5u(x, e 1 , . . . , e r ), 
Sg(x, e 1 , . . . , e r ), d5g(x, e 1 , . . . , e r ) and 6y(y, e 1 , . . . , e r ), 8u(y, e 1 , . . . , e r ), 

5g(y, e 1 , . . . , e r ), d5g(y, e 1 , . . . , e r ) are generally different, hence the func- 
tion forms of J£ and J* as functions of field quantities, their derivatives, and 
coordinates, are different. This does not ruin the covariance of our presenta- 
tion which is assured by the coordinate independence of eqn.(38). 

We have proven these for the case of GR.. However it is easy to show 
they are true for the cases of SR and pre-relativity dynamics. 
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According to Einstein, "[w]hat we call physics comprises that group of 
natural sciences which base their concepts on measurements; and whose con- 
cepts and propositions lend themselves to mathematical formulation." The 
aim of all natural sciences is to search for the objective laws of nature. 
In order to measure physical quantities or to describe physical processes in 
terms of mathematics, one needs first to choose a reference coordinate sys- 
tem which depends on the observer's subjective will. Therefore, the ways of 
performing measurements and the formulations of physical laws, expressed 
in all reference coordinate systems should take the same form. This idea 
has guided Einstein from pre-relativity physics in Newtonian absolute space- 
time, to SR in Minkowski spacetime, and finally to GR in curved space- 
time. In particular, to get the proper motion equation one has to use an 
invariant Lagrangian C(g(x), dg(x), d 2 g(x),u(x), 5w(x)),such that the action 
A[g,u] = J x(n) d 4 xy/-\g(x)\£(g(x), dg(x), d 2 g(x),u(x), du(x)) is an invari- 
ant under general transformation of coordinates. The least action principle 
would lead to a motion equation independent of coordinate systems. The gen- 
eral covariance of a physical theory means the operations of measurement, 
and the mathematical expressions of physical laws should be independent 
of reference coordinate systems. It does not mean that functions of coordi- 
nates, field quantities, and their derivatives, in different coordinate systems 
but with the same function form, must represent the same physical geomet- 
rical object, as is generally accepted. This is not always true as has been 
shown above. So, Noether's conserved current (38) 

J K [u(x), du(x), g(x), dg(x), d 2 g(x); 5x,5u(x), Sg(x), d5g(x)\ =: J£ and 

J K [u(y),du(y),g(y),dg(y),d 2 g(y); 5y,Su(y),5g(y),dSg(y)] =: J* are com- 
ponents of the same vector field on spacetime; while the canonical energy- 
momentum (41) Tp(x) = r K p (u(x), du(x), g(x), dg(x), d 2 g(x)) and r K p (y) = 
Tp(u(y),du(y),g(y),dg(y),d 2 g(y)) are not components of the same tensor 
field on spacetime. The long standing pseudotensor problems like (61), are 
only a result of misreading Noether's theorem; and the non-locality difficulty 
of gravitational energy-momentum does not really exist at all, had we read 
Noether's theorem properly. They are results of mistaking different physical 
geometrical objects as one and the same. 

It is important to distinguish a general law of nature and a concrete phys- 
ical object or process. The former is independent of coordinates, while the 
latter might looks different to different observers (from different coordinate 
systems). Yet we still can observe and describe a concrete instance in a way 
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independent of coordinate systems. 

Now we are in a position to address the problems of conservation of 
energy-momentum and the gravitational energy-momentum in GR. These 
will be done in a later paper. 

Acknowledgement 6 / am grateful to Prof. Zhanyue Zhao, Prof. Shihao 
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A Appendix 

The Lagragian density of classical fields L, is a function of the coordinates, 
field quatities, and their derivatives of up to the iV-th order. However, be- 
cause not all the variables are independent, such as d^d v u^(x) = dvd^u^x) , 
9a/3{ x ) — 9/3a( x ), etc., there are infinitely many different function forms for 
L. This causes indefiniteness of derivatives, such as — ® e , - L, » 9 , If 

we drop all the redundent variables, then the Einstein summation convention 
can no longer be used, and the expressions will become awfully complicated, 
especially for a large N. In order to keep the formulaes neat, physists usually 
treat it in a different way. Here we will illustrate their method by using the 
lagrangian density for vacuum Einstein's equation, R (Ricci's scalar curva- 
ture). 

R is a function of 16 g a p(x)% 64 d^g al3 (xys, and 256 d tl d u g al3 (xys. Be- 
cause g a p{x) = g/3a(x), and d^d u g a/3 (x) = d^^g^x), there are only 150 in- 
dependent variables among them. We will choose 10 g Q( g(a;)'s, 40 d^ga^xYs, 
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and 100 d /x d l/ g a ^(x) , s (a < (3,fj, < u), for the independent variables. As a 
function of 336 variables (As a function defined on a 336-demensional do- 
main), R can take infinite different forms, say, <p,ijj, . . .When restricted to 
the 150-dimensional "sub-domain" D, all of them are the same function of 
150 variables. 

R\d = ¥\d = i/)\d = ••• (Al) 

Substituting the 150 independent variables for all the variables in (p,ip, . . ., 
we get a unique function 

Ri,g a pi,x),d^g a p(x),d^d v g a p(x)), (a < (3 , fi < v) (A2) 

Substituting \(g a p(x) + g Pa (x)), \(d^g aP (x) + d^g^x)), and \{d^d v g aP {x) + 
d^d v g Pa (x) + d v dng a p{x) + d u d^g Pa (x)) for g afS (x), d^g a/3 (x), and d^d u g aP (x) 
in R, respectively, we get a unique function of all 336 variables, denoted 
by R(g(x),dg(x),d 2 g(x)). This "standard" R(g(x),dg(x),d 2 g(x)) has the 
following property. 

dR OR 1 dR 



dg a p{x) dg Pa (x) 2dg a/3 (x) 
dR dR 



, a < 



dg aa (x) dg aa (x) 

OR OR 1 dR 



(A3) 



dd^gapi 


» 


dR 




ddf,g aa { 


x) 


dR 




dd^d v g a 




dR 





dd^g Pa (x) 2dd lx g a/3 (x) 

dR 
dd^g aa {x) 

dR dR 



, a < (5 



(A4) 



ddpdv9p a (x) dd u d^g a p(x) 
1 dR 

a < f3,/i<is (A5) 



ddud^gp^x) Add^gapiyx) 
dR dR 1 dR 



dd^gaaix) dd^gaaix) 2 dd lt d v g aa {x) 
dR dR 1 dR 



dd^d^g aP (x) dd^d^g aP (x) 2dd lx d^ l g a p{x) 

dR dR 
dd^g aa (x) ddpdngaaix) 
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,H<v 
,a < (3 



(A6) 



When calculating the derivatives of R, we pretend that all its 336 variables 
are independent. Thus the indefiniteness problem no longer exists. 
From (Al), we have 

Slp\ d = 6ip\ D 

While 

dip dip dip 

r v-^ dip >r-^ d(f dip . . , . 

/x / dip 



a dd^gaaix) ^a</3 dd lx g a p{x) 

— 7 ; )Sd»gap(x) + V — ^ r Sd^gaaix) 

dd^gp^x) ^a,n dd^gaaix) 

dip dip 
^a<fs,^ dd^g a p(x) dd^d^gp a {x) 
dip dip 
+ ^ a ,ix<^dd^d u g aa (x) + dd u d^ 1 g aa ( y x)' >Sd ' ldu9aa(yX " > 

+ y ( *t + *e. + 

^ a <a^<v dd^g^x) ddvdpgapix) 

+ , , )5d,d u g aP {x)]\ D (A7) 



dd,j,d v gp a {x) dd u d^gp a (x) 

Because all the variations on the RHS of (A7) are independent, we get 

( dy? + dip = dj> + djj 
dg a p(x) dgp a (x) ° dg a ^(x) dgp a (x) ° 

( d( p + ft? )i = ( w + gV; )i 

dd^g^x) dd^g Pa (x) D dd^g^x) dd^g^x) D 

dd^g^x) ddvd^gapix) dd^g^x) dd v d^g Pa (x) D 
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^ dd p d u g aj3 (x) dd v d p g aj3 (x) dd^g^x) dd^g^x) ° 



This tells us, say, 



dg aj3 (x) 11 dg aj3 (x) dg a/3 (x) 

-5d^g a ^{x))\ D = ( — — ^ ^ a<3 (x))| D 



^dd p g aP (x) dd p g a/3 (x) 

#y^,0 a/ 3(x))|£> = ( a» ^ 7 ^ dpdvgapix)) 



^dd^d v g aP {x) ddnd v g a p(x 

dtp 



8dpd v g aP {x))\ D (A9) 



and even more, such as 



dd K g aP (x 
dip 



-){8g aP {x) - 5x p d p g aP {x))\ D 
■)(8g a p{x) - 5x p d p g a p(x))\ D 



dd K g a/3 (x) 

where R(g(x),dg(x),d 2 g(x)) is the "standard" expression for R, and 
ip(g(x),dg(x),d 2 g(x)) is any expression from (Al). 



B Appendix 

Proposition 7 Lei 

j« =■ ( dR - d dR - r (x) dR u 

x ' [ dd K g a p(x) "dd K d,g aP (x) ^ ) dd K d p g^(x) ) 

dR 

(Sg a p(x) - 8x p d p g a p(x)) + — - - d^Sg^x) - 5x p d p g al3 (x)). (Bl) 



Then 

Ly ~ dx^ x 



= (B2) 
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Proof. 



3 11 

-^d a g^dpg ar) - -d^g a/3 d v g pa - -d p g a( d v gf3 a - d a gp p d^g ar) ) (B3) 
f)R 1 

" ! = ^9^ + 9^9^) (B4) 



dd K d^g al3 2 

Note that g 9 g fl g is a (4, 0)-tensor, symmetrical for and for (a,/3). 

f ' = [g aK g^ + g^g<g^ + - <7<V V" 



-\g<*P g Wgt,v _ g< - g aR g^g m - g a ^g pK g m ] 
= :d lx g^B Ka ^ (B5) 



1 

OOKOfj.ga/3 Z 

+^(g aK d^ + g^d^) - g a ^ 



+d,ga-\9 aK 9^9 m - \<f*9 K <r + 9^9^9 m ] 
= : d^A"*™" (B6) 

Note that A Ka ^ and B Ka ^v are (6, 0)-tensors.Let 



= -\g a ^x)g^{x)g^{x) + ±<r(x)gK(x)gr(x) 
-\sr*(x)tf*{x)gr{x) - \g m {x)s^{x)sr{x) 

+ \g a Kx)9^{x)g^{x) + g aK (x)g^(x)g^(x) - \g af) {x)g^{x)g^{x) (B7) 
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Then 



0*) = t^ T (,) + A(gl) (B8) 



d^g^x) 



i OR d dx^ dR 

~ K '^ y) dd^g a ,^y) + ^W^O^g^iy/ 9 ^ 

dx K dR dx a dx 13 d dx a dx 13 - 

— OR — Ox K 

dx K r „i dx K 

= ^r^K + ^^ rest 
ay dy K 

d . dx^ dR - 

dx» dy» ddn'd^g^piy) 

We are going to show that rest vanishes. Its first term is 

= [g a ' e (y)g"' K ' (y)g" w (y) - \g a ' e (y)g" w (y)g K '"' (y) 
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where 



-g a ' K/ (y)g^'(y)g^'(y) + \g a ' p '{y)g^\y)g^\y)] * 

d 2 x^ dy^ 



gt"v>(y) s ga'p'(y) 



= g a ''\y)g"' K '{y)^{^)~8gMy) (A) 

1 f) 2 T^ Bill 3 ' - 



-g^\y)g^\y)^{^gMv) (A) 
+\^(v)g""(v)^(^)hMv) ((C) (bio) 



The second term is 



= (y)/ v (y)/"' (y) - ^ (y)/"' (y) <?" v (y) 



-^(v^W^i/) + ^ Q ' /3 '(y)/ /1 '(y)/' ? '(y)] x 



dy^'dy v ' dx v 



gev"(y) s 9a'^'(y) 



g"' K \y)g^\y)Q^%fhMy) (d) 
-^^'(y)^'(y)^^^W(y) (B) 



-y a ' K/ (y)/ v (y)^(^)^ r (y) (A) 
+^ a,/3 '(y)y^'(y)^(^)W(2/) (c) (bh) 
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The third term is 



d dx^ dR - 

= ^(^V^fe) (A) 



The forth term is 



-^(^7)^^ K ^^(y) (C) (B12) 
dR d 2 x a dy al - 



dy^'dy 01 " dx a 



a^l^V'Wfe) (B) 



s^i2l fyOL .Oil 



The last term is 



Sg a// 3/(y) 



dd K >d li >g a 'p"(y)dyi i 'dyP a dx? 

' V"V v - ^v v )W? 



dy^'dyP" dx^ 

'"'V^'W?/) (D) 



ch/' 3 " cte' 3 ' 

+ ww^ g ^ 9 SgMv) (B) (B14) 

All the terms marked (A) cancel each other, all the terms marked (B) cancel 
each other, etc. Therfore we get 

4 K = (B15) 

That is (B2). ■ 
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Abstract 

Based on a general variational principle, Noether's theorem is re- 
visited. It is shown that the so called pseudotensor problem of the 
gravitational energy-momentum is a result of mis-reading Noether's 
theorem; and in fact, all the Noether's conserved quantities in general 
relativity are scalars. It is also shown that the non-localizability of 
gravitational energy-momentum can not be attributed to the equiva- 
lence principle, by using a counter-example. As a direct consequence 
of variational princle, a generalized Hamilton-Jacobi equation for the 
Hamilton's principal functional is obtained. 



1 Introduction 



Ever since Einstein introduced the gravitational energy-momentum complex 
to keep the law of conservation of energy-momentum alive in general relativ- 
ity (GR) in 1918, the so-called pseudotensor problem of various gravitational 
energy-momentum complexes has been troubling relativitists [1]. One of the 
consequences of pseudotensor property is the non-localizability of gravita- 
tional energy-momentum. It is attributed to the equivalence principle. Some 
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relativitists think it is inherent in the theory of GR. However, not all relativi- 
tists accept this viewpoint. It is widely accepted that the proof of the posi- 
tivity of the total gravitational energy both in spacial and null infinity is one 
of the greatest achievements in classical GR in the last 30 years. This success 
inspired efforts to find quasi-local gravitational energy-momentum. However, 
finding quasi-local gravitational quantities has proven to be surprisingly dif- 
ficult [2] . Experimentally, the gravitational wave carrying energy- momentum 
has not been tested directly [3]. The aim of this article is to explore the pseu- 
dotensor problem of the gravitational energy-momentum, in particular, the 
difficulty of non-localizability of the gravitational energy-momentum. The 
whole argument is based on the variational principle of dynamics. It is shown 
that the non-localizability of gravitational energy-momentum can not be at- 
tributed to the equivalence principle, by using a counter-example. It is also 
shown that the so-called pseudotensor problem of the gravitational energy- 
momentum is only a result of mis-reading Noether's theorem, and mistaking 
different geometric physical objects as one and the same; hence the non- 
localizability difficulty does not exist at all. In fact, all Noether's conserved 
quantities in GR are scalars. As a direct consequence of variational prin- 
ciple, a generalized Hamilton- Jacobi equation for the Hamilton's principal 
functional is obtained; and it will be used for the quantization. of fields in 
GR in a later work. 

The history of modern physics has proven that the variational principle 
approach to dynamics is not only an alternative and equivalent version to the 
naive, intuitive approach, but also yields deeper insights into the underlying 
physics. For instance, it is hard to imagine that the statistical mechanics 
could have been established without using the concept of phase space, and 
that the quantum mechanics could have been established without using the 
concept of Hamiltonian. Therefore, we will base our argument on a gen- 
eral variational principle for the classical field. It might be for the similar 
consideration that soon after Einstein proposed his general theory of rela- 
tivity, Hilbert made the first succesful attempt to get Einstein's equation 
by using the least action principle. The Lagrangian being used for vacuum 
Einstein's equation, (16irG)~ 1 R, is the only independent scalar constructed 
in terms of the metric field and its derivatives of no higher than the second 
order. However, because the Ricci scalar curvature R contains the second 
order derivatives of the metric field g^ v {x), which is now the dynamic vari- 
able along with the matter field, the least action priciple for Lagrangians 
containing only the field quantity and its first order derivatives does not lead 
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to Einstein's field equation. The generally accepted solution to this diffi- 
culty is adding the Gibbons-Hawking boundary term to the Hilbert action 
and keeping the least action principle unchanged[4]. But there is another 
solution to this difficulty, which was initiated by Hilbert and will be adopted 
in the present paper. The least action principle will be restated and the 
Hilbert action will still be used for the vacuum Einstein's equation. In order 
to show this is proper and natural, we will consider the variational principle 
for classical fields in 2)-dimensional spacetime with a Lagrangian con- 
taining the field quantity and its derivatives of up to the iV(^ l)-th order. 
In our opinion, non-local interaction is not acceptable, so we assume that the 
Lagrangian does not contain the integration of the field quantity. 

There have been varied versions of variational principle and Noether's 
theorem in the literature [5], and different notations have been used by dif- 
ferent authors. For the readers' convenience, and for the consistency of the 
reasoning, we start with presenting a general variational principle for classi- 
cal fields in 2)-dimensional spacetime with a Lagrangian containing the 
field quantity and its derivatives of up to the l)-th order. In section 2, 
the Lagrangian formalism is presented, and a generalized Hamilton- Jacobi 
equation for Hamilton's principal functional is obtained. In section 3, the 
Noether's theorem is rederived. Then, in section 4, these general results are 
applied to the specific case of general relativity, especially Noether's theo- 
rem is applied to get quite a few conservation laws in GR. It is noted that 
most of the conserved quantities obtained here are not tensors like various 
gravitational energy-momentum complexes. By using a counter-example, we 
show that the non-localizability of gravitational energy-momentum can not 
be simply attributed to the equivalence principle. Then Noether's theorem is 
revisited in section 5, and it is shown that the so-called pseudotensor problem 
of varied gravitational energy-momentum complexes is a result of taking dif- 
ferent geometrical physical objects as one and the same; as a matter of fact, 
all the Noether's conserved quantities in GR are scalars; and the non-locality 
difficulty does not really exist. 

These results will be used to explore the energy-momentum conservation 
and the gravitational energy-momentum in GR in a later paper. 
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2 A general variational principle for classical 
fields 



2.1 Lagrangian formulation 

First, we present a useful mathematical formula for the variational principle, 
which does not rely on physics. 

Suppose {$b : R n — ► R (or C) | B = 1, 2, . . . , /} (n ^ 2) are smooth 
functions, and function L = L(x, $(x), d<&(x), d 2 &(x), . . . , d N <&(x)) is smooth 
with respect to all its arguments. It is easy to show just by using Leibniz's 
rule that (the Einstein convention is used for coordinate indicesjj 

f - 81 

f N f N-l 

B=l X=0 Ai A X _«v ; B=l z=o 

5d Vl ---d Vz ® B {x) (-!) r <V--3 



v _ u d<9 A d Ml ■ • • d^du, ■ ■ ■ d Uz $ B (x) 

f V dL 



/ AT-1 

+ ^ E E ^ " " " d» z ®B{x)K BXv ^{x) (1) 



1 In order to avoid the indefiniteness of derivatives such as gg 1 g 2(S , B L, and in order to 
keep the formulaes neat, it is assumed in the present paper without loss of generality that 

L{x, ${x),d$(x), d 2 1>(x), d N $(x)) 
= : L(x a ,$ B (x),d ai $ B (x),d ai d a2 $ B (x),. . . ,d ai d a2 ■ ■ ■ d aN $ B (x)) 

= L(x a , $ B (x),d a $ B (x), - Y, PeS2 d ap{1) 9 ap{2) ^ B (x), . . . , 



We will pretend that all the cross derivatives are independent variables of L, when calcu- 
li 

dd 1 d 2 'S>E 



lating »„ i . L, etc. See Appendix A for the details. 
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when $ B (x) 1 — > $b(x) = ®b(x) + 5$ B (x),VB = l,...,f. Consider the 
functional F of the following form 

= I d n xL{x,®{x),d${x),d 2 ®{x),...,d N <S>{x)), (2) 
Jn 

where f2 is an open subset with a compact closure of 1R™. When the arguments 
Q B (x) \/B = 1, . . . , / change slightly, the variation of functional F is 



5F[*] = [ d n xJ2J25d Xl -..d Xx ^ B (x)—- 



d Xx $ B (x) 



, f N-l 

+ / d n xd x j2J2 Sd ^--- d ^ B ^ KBXui " uz ^ 



B=l Z=0 



, f N-l 

«/ o 1 v n 



B=l Z=0 



This can be easily obtained by using eqn.(l) and the Stokes theorem. The 
derivatives of functional (2) is defined as follows. 



Definition 1 If the change of the functional (2) can be expressed as 

r f 

F[$ + 6<f>] - F[$] = / <rx^26$ B (x)D B [$,x\+o[6$\ 



(4) 



where D B [&,x] is a functional of $ varying with x, and o[5§] is a higher 
order infinitesimal of 5$, when | B — 1, 2, . . . , /} change slightly while 
the boundary values of $, cM>, . . . , d N ~ l <& are kept fixed, then F is called dif- 
ferentiate at $ 7 and D B [Q,x] is called the derivative of functional F with 
respect to at $ and point x, and denoted by 
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»™ -a 

Let us now apply the general formula (3) to the action functional of 
classical field {$# : M — > R (or C) | B = 1, 2, . . . , /} 

A[$]= f d n xL(x,^(x),d^(x),d 2 ^(x),...,d N ^(x)) (6) 

Jx(Q) 

where Q is an open subset with a compact closure of the spacetime mani- 
fold M, x(Q) C M. n is the image of Q C M under the coordinate mapping 
x : M— > R n and L is the Lagrangian of the field. We get the difference be- 
tween the action functionals over Q of two conceivable movements (two paths 
allowed by the constraints) close to each other 



r f N fit 
6A[*]= / (Ti^^EH^-^as d~^lx) 

Jx(si) bT-l £r oo Xl - ■ ■ o Xx ^b{x) 

r f N-1 

+ / d n xd x Y,Y, 8d "i'~ dv z ®B(x)K B ^-»*(x) 
Jx(n) B=1 Z=Q 

r f N fit 

, f N-1 

+ / ds x {x)Y J Y. 5d ^--- d -^B{x)K BX ^- Vz {x) (7) 
Jx(an) B=1 Z=Q 

Equation (7) suggests that for all AT ^ 1, n ^ 2 the least action 
principle read as follows. 

For any spacetime region Q, among all possible movements in Q with the 
same boundary condition 

8<£>\ dQ = 0, 5d$\ dn = 0, . . . , Sd^^U = 0, (8) 

the real movement (the path allowed by physical laws) corresponds to the 
stationary value of the action over Q. 

Combining eqns.(7), (8), one obtains the field equation (Euler-Lagrange 
equation) satisfied by the real movements 
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This Lagrangian formalism applies to Newtonian dynamics, dynamics in 
special relativity (SR) and dynamics in general relativity. Specific covariance 
is the heritage from the Lagrangian being used. It is not the heritage from 
this general formalism. 



2.2 Hamilton's principal functional and Hamilton-Jacobi's 
equation 

Let us consider the difference between actions over spacetime region Q of two 
real movements close to each other. Using eqns.(7) and (9), one gets, for real 
movements 

, / N-l 

5A[<S>] = / ds x (x) 5d »i ' • -du z ^B(x)K B ^- Vz (x) (10) 

Jx(dn) B=1 z=0 

From eqn.(10), one sees that the action over a spacetime region Q of a 
real movement is determined by the closed hyper-surface dQ, and 
d$\gn, . . . , 5 Ar_1 $|an- It will be called the generalized Hamilton's princi- 
pal functional and denoted by 

S = S[dQ, Siangan, ... , (11) 
Re- write eqn.(10) as 

, / N-l 

SS= ds x {x)Y J Y, 6d ^--- d -^B{x)K BX ^-^{x) (12) 

Jx(dn) B=1 z=0 

Note that when ^slan is given, only one of the n derivatives d\$s\dn 
(A = 0, 1, . . . , n — 1) is independent; when dx&sldn is given, only one of 
the n derivatives d tl d x ^B\do. (/-* — 0, 1, . . . , n — 1) is independent; and so on. 
Thus for a given suffix B, only N items ( ^s\dn, ° ne of dx^slan, ■ ■ ■ , ° ne of 
d\i • • • d\ N -^B\dn {Xj = 0,1,..., n-l) ) are independent. 

In order to formulate the generalized Hamilton-Jacobi's equation, one 
needs a new type of functional derivative [6]. 
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Definition 2 Let E be a closed hypersurface in W\ ^ a function defined on 
M. n , and F = F[E,\I/|e] a functional o/E and The functional derivatives 
are defined as follows. If the variation of F can be written as 

6F\E, = J ds x (x) {y[£, x^E^x) + Z[E, x] A o^(x)} (13) 

where Y[H, x]* and x] A are junctionals of E and \l/|s, which 

vary with x, when E i — >■ E, and ^ i — >■ \l/ 7 taen £/iey are called the functional 
derivative of F with respect to E M (x) and ^(rc), and denoted by 



(5F 



respectively. 

Hence we have 



5F[E,tf| s ] = ^ds A (x) 



<yE"(x) + 



The hypersurface E is given by the parameter equation 
The £E M (:r) and S^f(x) in eqn.(13) are respectively 



6V(x) 



(14) 



(15) 



(16) 



5W(x) = W(9 l ,9 2 , 9 3 ) — E M (^ 1 , 9 2 , 9 3 ), 
6V(x) = #(E"(0\0 2 ,0 3 )) -^f(W(8 1 ,9 2 ,9 3 )). (17) 

Now, from eqns.(12) and (15) we get [dQ in (12) is E in (15)] 



5S 



Sd ui ■ ■■d Vz $ B (x 



}j = k BXvi - Vz (x), 



0,...,N-l,B = l,...,f. 



(18) 



Following the evolution of one real movement and observing the change of 
its action with E, 
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6S= f ds x {x)L5 x a 5T 1 (7 {x) 

Jx(dQ) 

f N-l , \ A 

+ ^■■■W) (19) 

B=0 Z=0 v 1 z v y y 

We get the generalized Hamilton- Jacobi's equation. 

A / N-l 

+ [L<£ - £ J] d ^ ■ ' ' (*)] = (20) 

B=0 z=o 

It will be shown in the next section that the expression in the bracket is 
just the canonical energy- momentum t*(x), hence we have 

For a dynamic system with finite degrees of freedom, we have the Hamilton's 
principal function S(q 1 , . . . ,q h ,t), and the Hamilton- Jacobi's equation is a 
partial differential equation || + H(q,^,t) = 0. While for the classical 
field, we have the Hamilton's principal functional, and the Hamilton- Jacobi's 
equation is a set of differential equations for functional. It plays an important 
role in canonical quantization for fields. 




3 Noether's theorem 
3.1 Re-deriving the theorem 

There have been varied versions of Noether's theorem, and different notations 
have been used. To avoid confusion, here we redrive Noether's theorem for 
classical fields in n-dimensional spacetime with a Lagrangian containing field 
quantity and its derivatives of up to the iV-th order, in terms of coordinate 
language. 
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Theorem 3 If the action of classical fields over every spacetime region Q 
with a compact closure 

A[$]= [ d n xL(x,<5>(x),d<5>(x),d 2 <5>(x),...,d N $(x)) 

Jx(Q) 

remains unchanged under the following r— parameter family of infinitesimal 
transformations of coordinates and field quantities 

x x i — > x x = x x + Sx x , 
5x x =: 5x x (x, e 1 , . . . , e r ), \e l \ < 1, 5x x (x, 0, . . . , 0) = 

$ B (x) 1 — > $b(x) = $b(x) + 5<Z> B (x), 
5$ B {x) =: 5$ B (x, e 1 , . . . , e r ), \e l \ < 1, 6$ B {x, 0, . . . , 0) = (22) 

£/ien there exist r conservation laws. 

Proof. The small change of field quantity 5&b{x) can be devided into two 
parts, the part due to the small change of its function form and the part due 
to the small change of its arguments respectively. 

5<5> B (x) = 5<5> B (x) + 5x a d a $ B (x) (23) 

Similarily the small change of derivatives of field quantity 5[d\ 1 d\ 2 ■ ■ ■ 8 Xx ^b(x)] 
can be written as 

S[d Xl dx 2 ■ ■ ■ d Xx $ B (x)} = d Xl dx 2 ■ ■ ■ d Xx M B (x) 

+5x°d a d Xl d X2 ---d Xx $ B (x) (24) 

The variation of the action can be devided into two parts. One is due to the 
small change of the integration domain x(Q) i — > x(Q) in M n , and the other 
is due to the small change of the integrand 

L(x, $(x), d$(x), a 2 $(x), . . . , d N $(x)) 
i— >■ L(x, d$(x), d 2 $(x), . . . , d N $(x)) 
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6A[$] = [ ds x (x)5x x L 
Jx(dn) 

f n 



r bt 

+ / ^EE 99 9-^) d ^ ■ ■ ■ d >J*B{x) 

Jx(Q) B =1X=0 ^ ■■■V^x^BK.X) 

f N ~j 

B=l X=0 Ml °>x VB W 

, / N-l 

+ / ^^[^L + ^^a^-..^^^)^^—^)] (25) 
Jx(n) B=1 z=0 



L 



The first integrand at rhs vanishes for real movement, hence the second 
integral does too. One gets the following continuity equation due to the 
arbitrariness of Q. 

f N-l 

d x [5x°5 x L + £5>i-^ (6$ B (x) - 5x°d a <f> B (x)) K BXvi - Vz (x)\ = 

B=l Z=0 



J 

Jx( 



or 

/ N-l 

ds x (x){5x°5 x L + Y,Yl 
ix(dn) B=1 z=0 

[d Vl • • • 3„ z (6$ B (x) - 5x°d a <$> B {x)) K BX ^- Vz {x)]} = (26) 

Noting that both 5x a and 5§ B (x) depend on r real parameters, one gets from 
eqn.(26) r conservation laws. 

/ JV-l 



B=l Z=0 ^ ' 

+ ^| e=0 ^L} = 0, Va = l,...,r (27) 



3.2 Conservation law due to "coordinate shift" invari- 
ance 

In this subsection, we restrict our discussion to those cases, such that the 
Lagrangian does not manifestly contain coordinates and the action of classical 
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field over any spacetime region 

A[$]= f d n xL(^(x),d^(x),d 2 ^(x),...,d N ^(x)) (28) 

Jx(Q) 

remains unchanged under the following "coordinate shift". 

5x a = e\ 8® B (x) = 0. (29) 

And eqn.(27) reads 

d x r x a (x) = 0, (30) 

where 

/ N-l 

t*(x) =*!l-tE 9 A • ~%>**b{x) K b ^—{x) (31) 

B=l Z=0 

is usually called canonical energy-momentum. 

Comparing eqn.(20) and (31), we get eqn.(21), the generalized Hamilton- 
Jacobi equation. 

The formalism presented so far is good for any classical field in n-dimensional 
spacetime with a Lagrangian containing field quantities and their derivatives 
of up to the iV-th order, no matter it is Galileo covariant, Lorentz covari- 
ant, general covariant or without any covariance. Therefore, we can not say 
eqn.(21) is covariant even though it looks like so. The covariance of the for- 
malism is the heritage from the Lagrangian being used, not the heritage from 
the formalism itself. 

4 Variational principle approach to general 
relativity 

Let us apply the general results obtained above to the classic fields in GR. 
We will consider the case where the matter field is a (1, l)-tensor field. The 
results can be readily generalized to any (r, s)— tensor matter field. For the 
dynamic system, (1, l)-tensor matter field u|(x) plus the metric field g a p{x), 
the Lagrangian and the action over any spacetime region Q are 
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L(g(x),dg(x),d 2 g(x),u(x),du(x)) 

y/-\g(x)\[C(g(x),u(x),Vu(x)) + -J—^R] = L M + L G , (32) 



IQttG 

A[g,u] = [ d 4 x^-\g(x)\[C(g(x),u(x) : Vu(x)) 
Jx(n) 

+ T ^R} = A M [g,u} + A G [g] (33) 

where R is the Ricci scalar curvature, C(g(x), u(x), Vu(x)) is the Lagrangian 
for matter field obtained from the Lagrangian in special relativity 

£(r/, u(x), du(x)) by replacing the Minkowski metric 77 = diag(—l, 1, 1, 1) 
with g a/3 (x), and replacing the partial derivatives d^u^(x) with the covariant 
derivatives V^u°(x). 

4.1 Einstein's field equation 

The Euler-Lagrange equation, Eqn.(9) now reads 

^=«4b- v ^ ] =° (34) 

W) = - \ R ^ X) - 8 * Gr * 1 = ° (35) 

where T a/3 is the energy-momentum tensor of matter field, which is a sym- 
metrical (2,0)-tensor. 

rpafs = 2 j4M _ N , dC 

~ V=\M\ S 9<#(x) ~ 9 {X)L + 2 dg a ,(xy u ^ x) 

, dc avAttg(s) a£ sv^f (x) 



5Vam|(x) dg a p{x) ' v dV^u^x) dd x g a p(x) 
2<9> — — ^— — — - 2— — 



dV^u^x) dd x g aP (x) <9V dd x g a p(x) 
T al3 (u(x),du(x),g(x),dg(x)) (36) 
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All the pre-GR dynamics discuss how matter behaves in spacetime with a 
given metric, that is, Newton's absolute space-time or Minkowski space. All 
of them can not explain why there are inertial reference systems and non- 
inertial reference systems, and why the inertial mass equals the gravitational 
mass. G.R. is unique. It discusses how the matter motion determines the 
bending of the spacetime and how the matter behaves in a curved spacetime. 
Hence the dynamical variables are both the matter field u®(x) and the met- 
ric field of spacetime g a p(x). It means that they are both to be determined 
simultaneously by solving the equation of motion. It does not necessarily 
imply that the metric field of spacetime g a p{x) is matter (in its philosophic 
context, in the most general meaning of the word) like, say, the electromag- 
netic field, the spinor filed, etc. It does not necessarily imply that the metric 
field of spacetime g a p{x) carries energy- momentum. 

4.2 Noether's theorem for classical field in GR 

The Noether's conservation law, or the continuity equation (26) , now reads 



J K [u(x),du(x),g(x),dg(x),d 2 g(x);5x,Su(x),5g(x),dSg(x)]} = (37) 



d 

dx K 



{V-\9(x)\J K [u(x),du(x),g(x), dg(x), d 2 g(x); 5x, 5u(x),5g(x),dSg(x)] j 







or 




where 



J K [u(x),du(x),g(x),dg(x),d 2 g(x);Sx,Su(x),5g(x),d5g(x)] = 




g 



£ a (x)u°(x) 
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dC dC - 

1 AR^H^^-dl dR 



r U x ) l^ a ~ _ f J &9<#{x) + aaa „ ^ dpSgapix)} (38) 



167rG L " p dd K g a/3 (x) dd K d^g af} (x) 

dR - dR 
' dd^ga^x) al3 dd K d li g a p(x) 

It is worth noting that the form of function J K is independent of coordinate 
systems, and the arguments of function J K are u(x), du(x), g(x), dg(x), d 2 g(x); 
and 8x, Su(x), Sg(x), d5g(x). 

4.3 Conservation laws in GR 

4.3.1 Conservation law due to "coordinate shift" invariance 

Action (33) remains unchanged under the following "coordinate shifts". 

5x p = e p , 6ul(x) = 0, 8g aP (x) = 0. (39) 
In this case, eqn.(26) reads 



dxly/lgWKW] = o, (40) 

where 

T p( x ) = i~p(u(x),du(x),g(x),dg(x),d 2 g(x)) 

dC J7TM^ - ^ + V^Juy a {x)4{x) 



dV K u\{x) p ^ ' p TdV K ul(x) 
dC dC 



av K ^(x) y v ' ^ J dVpuiix)- 

dC I 

+ aw a ( , 9 iK {x)u e Ax)}d p g a p{x) 



dR dR . . dR . . . 



'dd K g a/3 (x) M dd K d^g a/3 (x) vp dd K d li g a p(x 
dR 
dd K d^g a/3 (x 



+ ^ fl : .t^ dpg^x) - itfj (41) 
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is usually called canonical energy-momentum tensor. It is worth noting that 
the arguments of function Tp(u(x),du(x),g(x),dg(x),d 2 g(x)) are all field 
quantities and their derivatives. Notice that "coordinate shift" eqn.(39) 
is not an invariant concept under general coordinate transformation. This 
explains why t x (x) is not a tensor under general coordinate transformation. 
We will get back to this problem later. 



4.3.2 Conservation law due to "4-dimensional rotation" invari- 
ance 

The action (33) remains unchanged under infinitesimal "4-dimensional ro- 
tations" (Lorentz transformtions), which form a 6-parameter family of in- 
finitesimal symmetry transformations 



x m _> W = L p x u , Ut = 8$ + A£, \K\ « 1, V ^A X = A M „, = -A„„, 
^n xp x a -n Xa x p )A„„. 



5x x = A*a? = hrj , *x ff -rfx^Apr, (42) 



8u\{x) = Aj«?(x) - iqv°(x) 



1 



= A pa -[ v ep u°(x) - rf°u'{x) - 5\rf»>u e v {x) + S^ufa)], (43) 

Sg a p(x) = ApJ^-Slri^g^x) 

+ SirTg^x) - 5X P 9au{x) + 5 p ^g au {x)} (44) 
In this case eqn.(26) reads 

d 

dx K 

where 



{yTW)\M^} = (45) 



2M Kpa ( y u( y x), du(x), g(x), dg(x), d 2 g{x)) 
= C5 K x (rj Xp x° - r^af) + dVKU e {x) {lv 9p ^) ~ ^(x) 
- 5^ p u d v (x) + S^ufc)] - {r J Xp x (T - r] Xa x p )d x ul{x)} 
+ ~ h \, , g 6a (x)u P Ax) + - ftl - g ea (x)u K Ax) 
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dC 



dC 



: <f tt {x)u\{x) H 



dC 



g" a {x)u\{x) 

g^(x)4(x)]x 



dVpuHxY ' dV a u^{x) 

{[-^g^{x] + SZrTgtfix) - SfrTg^x) + 5 p V ™g a „(x)} 



- (v Xp x a ~ r} Xa x?)d x g a p(x)} + ^±^{RS K x (r^x ff - r^x") 

+ a ^ _V (x) ^ ) 

dd K g a p(x) 11 dd K d^g aP (x) ^ dd K d^g a p(x) 

([-*° a rrgtf(x) + SfrrgtfW - V P 9a,(x) + 8farg m {x)\ 

BE 

- - ^) w*)) + ^^^([-^^(-) 

+ KvTg^x) - SX P 9au(x) + 5 p ^g av {x)] 

-{r, Xp x° -v Xr7 x p )d x9aP (x))} (46) 

4.3.3 Conservation law due to "4-dimensional pure deformation" 
invariance 

The action (33) remains unchanged under infinitesimal "4-dimensional pure 
deformations", which form a 6-parameter family of infinitesimal symmetry 
transformations 

x m _> xf> = l£a7, l£ = % + AZ, |A£|«1, 
rj^At = A M „ = A^, (47) 

fe A = Ajx" = -{r] xp x a + ?7 A<T x p )A p(7 , (48) 
<5«f(s) = Aj«?(x) - AJ«5(x) 
= + i^«f(^) - ^V P <(*0 - ^<(x)]A pCT , (49) 

Sgap(x) = —K^g^x) 

+ W a g^) + V p g av {x) + S p V ^g au (x)]A pa (50) 
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In this case eqn.(26) reads 

d 



where 



dx K 



{y^(x)\N^}=0 (51) 



2N^(u(x),du(x),g(x), dg(x), d 2 g(x)) 

8C 
uv K u^yx) 

-Slrpfufa) - S^ulix)] - (r] Xp x a + r) X,T x p )d x u\{x)} 



1 f)r Fir 

c)C c)C 

lenK9 6 \x)4{x) - ^r-TTZK9^)4{x) 



p ua 



+ (v Xp % a + V Xa x p )dxg a p(x)} + ^{RS^x' - V x ° X o) 

_ d_R _ d dR ^ OR 

dd K g a p(x) P dd K d^g aP (x) V{1 dd K d^g a p{x) 

% 



+ 5 p ^g av {x) + (tf*x a + ^x p )d x g a p{x)\ 

+ 6 p V ^g au (x) + ( V xp x° + rj Xa x p )dxg a p(x)}} (52) 

4.3.4 Conservation law due to "scaling" invariance 

The action (33) remains unchanged under infinitesimal scaling transforma- 
tions, which form a 1-parameter family of infinitesimal symmetry transfor- 
mations 



18 



x x i — >x x = e e x x , |e| < 1, &r A = ex\ VA = 0,1,2,3 
<Jiif(a;) = 0, Sg a p(x) = -2eg a/3 (x) 
In this case eqn.(26) reads 



{v /r bMM = o 



(53) 
(54) 

(55) 



S\u(x)^x),g(x),dg{x),&g{x)) 



p dv K u\(xy dxp 



9£ n 9 a, . 

Tn-^x p — u e Ax) - - 



1 r <9£ g , . a. . 

5 {x)v%(x) 



2 l dW K u\{x) 



r)C r)C r)C 

fic fir 

<* a (x)u°(x) + ^ 0f J \x)uUx)}{2g a p{x) + x p d p g aP (x)} 



9V p u e K {x) 
1 



{R5> 



dR 



9 



OR 



KM- 



OR 



IQnG p 9d K g aP (x) ^ 9d K d fl g a p(x) 9d K d lx g af) (x) 



[2g a p{x) + x p d p g a/3 (x)] - 



OR 



99 K d fl g afS (x) 



9 fl [2g a/3 (x) + x p d p g afi (x)]} (56) 



4.3.5 Conservation law due to "skew-scaling" invariance 

The action (33) remains unchanged under infinitesimal "skew-scaling" trans- 
formations, which form a 1-parameter family of infinitesimal symmetry trans- 
formations 

Sx° = -eV, dx 1 = eV, 5x 2 = 0, Sx 3 = (57) 

" 1 - e 1 
9x] 1 + e 1 

TXT 7 ] ~ o oio 

1 



(58) 



Su° (x) = 0, 5u° l (x) = -2e 1 w°(x), 5u° 2 (x) = -e 1 ^), 5u° 3 (x) = -e^x) 
Sul(x) = 2e 1 ul(x), 5u\(x) = 0, 5u\(x) = t l u\{x), 8u\(x) = e^^x) 



Su J (x) = e u J (x), 8u\(x) = —e u\(x), 5u 3 k (x) = 0, Vj, k = 2,3 
$9oo(x) = Z^gooix), 6g Q1 (x) = 0, Sg 02 (x) = e 1 g 02 (x), 5g 03 (x) = e 1 g 03 (x) 
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Sg w (x) = 0, 6g n (x) = -2e 1 g u (x), 8g 12 (x) = -e 1 ^ 2 (x), Sg 13 (x) = -t l g n (x) 
Sgjo(x) = e 1 g j0 (x) : 5g jl (x) = -e 1 g j i(x), 5g jk (x) = 0, Vj, k = 2,3 (59) 
Substitute eqns.(57) and (59) into eqn.(26), we can get a conserved current 



Here we skip the expression of Jf (u(x), du(x),g(x), dg(x), d 2 g(x)). Similarily 
we can get conserved currents J% and Jg . 

4.4 Equivalence principle and non-localizability 

Because the symmetry group of classical field in GR is an infinite dimensional 
Lie group, we can get infinitely many conservation laws by using Noether's 
theorem. It is interesting to note that most of the conserved currents are 
not tensors. For instance, r p is not a (1,1) tensor, M Kpa is not a (3,0) tensor, 
N Kpa is not a (3,0) tensor, and S K is not a (1,0) tensor, etc. One of the con- 
sequences of the psedotensor property is the non-localizability. Relativitists 
attribute the non-localizability of the gravitational energy-momentum to the 
equivalence principle. The argument goes as follows. In a freely falling space- 
ship, the astronaut feels no gravitational field acting on him; while people on 
the earth see that the gravity is acting on him to keep him in the orbit around 
the earth. Hence the existence of the gravitational field at a spacetime point 
depends on the reference coordinate system adopted; non-localizability of 
gravitational energy-momentum is the direct consequence of the equivalence 
principle; and it is inherent in the theory of GR. This sounds very reasonable. 
But, let us examine the following example, Landau-Lifshitz's gravitational 
energy-momentum pseudotensor [7] 








(60) 



tr(x) 



[g» x (x)g™(x) - g^(x)g Xa (x)] + g^\x)g^(x)[T^(x)T^(x) + 
K^n p (x) - n p (x)Y{ a (x) - TUx)T; p (x)] 

+g 1J \x)g^{x)[V%{x)V^{x) + T%(x)T%(x) - T» p {x)T{ a {x) 

-rLW»] +g xa {x)g pp {x)[v%{x)vi p {x) -r^(x)r^ p (x)}}. 
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When the energy- momentum tensor of matter vanishes t a/3 (x) = 0, the Ein- 
stein field equation has the solution g a f3{x) = )j a ». The spacetime is the 
Minkowski space, and the coordinate system (x° , x 1 , x 2 , x 3 ) is a coordinate 
system of inertia. In this coordinate system, 

t^(x) = 0,V0 ^ fi,u ^ 3. 

Let us switch to coordinate system y 1 , y 2 , y 3 ) =: (t, r, 9, if), such that 



t = x°,r=y/(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 , 

9 = cos -1 — =, 
^{x 1 ) 2 + {x 2 ) 2 + {x 3 ) 2 



ip = tan — . 

x L 

In this spherical polar coordinates, 

Note that the coordinate system (x° 3 ) and the coordinate system 

(y°, y 1 , y 2 , y 3 ) =: (t, r, 9, <p), unlike the above mentioned spaceship and earth, 
are not in relative motion; they share a common reference system. Yet for 
any spacetime point p G M,we have ^{x) = 0;and t 00 (y) < 0. This counter- 
example shows that the non-localizability of gravitational energy-momentum 
in GR, is far beyond what can be simply attributed to the equivalence prin- 
ciple. 

One of the main tasks of the present work is to show that these pseu- 
dotensor problems are just results from misreading Noether's theorem. 

5 Noether's theorem revisited 

Let us start with the simplest example. The canonical energy-momentum 
Tp(u(x), du(x), g(x), dg(x), d 2 g(x)) in eqn.(41) is not a tensor 

T*(u(x),du(x),g(x), dg(x), d 2 g(x)) ^ 

T*(u(y),du(y),g(y),dg(y),d 2 g(y)) . (61) 



dx K dy a ^ Xi f A ;i ; ; 

dy x dxP 
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This can be understood as follows. The 4-parameter families of infinitesimal 
transformations of coordinates and field quantities eqn.(39) and 



Sy" = e», 8u\(y) = 0, 8g aP {y) = 0. (39') 

do not correspond to the same 4-parameter family of infinitesimal diffeomor- 
phisms of spacetime M, when 

Q2 tx 

3 p e M, and < a, (3, a < 3, such that - y n _ L ^ 0, (62) 

ox a oxP 

Therefore, generally speaking, the conservation laws (continuity equations) 
corresponding to them are not equivalent to each other, and the conserved 
quantities corresponding to them are not the same geometrical physical ob- 
ject. Inequality (61) is a comparison between components of two different 
geometrical physical objects which have been mistaken one and the same. 

In general, infinitesimal transformations of coordinates and field quan- 
tities expressed in different coordinate systems with the same form (like 
eqns.(39) and (39')) do not correspond to the same infinitesimal diffeomor- 
phism of spacetime M. Conversely, the same r-parameter family of infinitesi- 
mal diffeomorphisms of spacetime M is described differently in different coor- 
dinate systems. For instance, if a 4-parameter family of infinitesimal diffeo- 
morphisms of spacetime M is described in coordinate system (x°, x 1 , x 2 , x 3 ) 
by eqn.(39), then it would be described in coordinate system (y° , y 1 , y 2 , y 3 ) 
not by eqn.(39'), but by 

which generally is no longer "coordinate shift" . 

It is generally accepted that in GR, the same physical geometrical object 
is expressed in different coordinate systems by functions of field quantities, 
their derivatives, and coordinates, with the same function form. (For in- 
stance, T K p (u(x), du(x),g( y x),dg( y x),d 2 g(x)) of (41) and r K p (u(y), du(y), g(y), dg(y),d 2 g(y)) 
are generally considered the components of the geometrical physical ob- 
ject, canonical energy- momentum.) But this assertion is not always true. 
For instance, corresponding to the same r-parameter family of infinitesi- 
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mal diffeomorphisms of spacetime M, the Noether's conserved current (38) 
J K [u(x),du(x), g(x), dg(x), d 2 g(x); 8x,Su(x),Sg(x),dSg(x)] , and its counter- 
part in coordinate system (y°, y 1 , y 2 , y 3 ), 

J K [u(y) , du(y) ,g(y), dg{y) , d 2 g(y) ; 5y, 8u(y) , 8g(y) , dSg(y)] , are not func- 
tion of field quantities, their derivatives, and coordinates, with the same func- 
tion form. (Because J K manifestly contains 5x, Su(x), Sg(x), dSg(x), besides 
field guantities, their derivatives and coordinates, while 5x, 5u(x), Sg(x), d5g(x) 
and 5y,5u(y),5g(y),d5g(y) are different functions of field quantities, and 
their derivatives, and coordinates.) Yet they are the same physical geomet- 
rical object, as will be shown in the following. 

To prove this, let us show that continuity equations written in different 
coordinate systems, but corresponding to the same r-parameter family of 
infinitesimal diffeomorphisms of spacetime M , are equivalent to each other, 
even though the conserved currents written in different coordinate systems 
are different functions of matter field, metric field, their derivatives, and 
coordinates. 

Let {A(e 1 , . . . ,e r ) =: A : M -> M | |e*| < 1, VI < i < r} be an r-parameter 
family of infinitesimal diffeomorphisms of spacetime manifold M onto itself, 
and A(0, . . . , 0) the identity mapping of M. For any chart (U, ip) of M, de- 
noting by ( 3 ) the corresponding coordinate system, we have 

(x°, x\x 2 , x 3 ) ^T^T" x x = x x + Sx\x; e\ . . . , e r ) =: x x + 8x x , (63) 
where ir x : (x°, x 1 , x 2 , x 3 ) x x is the projection operator, 

5ul(x;e\...,e r ) = : 8u\{x) = ^^(x) - ^-^(x), 
8g a p{x] e 1 , . . . , e r ) = : 8g aP {x) = -^^g^(x) - ^^-g av {x) (64) 
and the continuity equation (37) which is re-written for simplicity as 

' iV^WjlW = (65) 



dx K 



Switching to coordinate system (y°, y 1 , y 2 , y 3 ), which corresponds to chart 
(V, ip), we have 

y\ y\ y") ^ W y x = y x + 8y\y, e\ . . . , <T) =: y x + 8y, (63') 
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^(y; e 1 , . . . , e r ) = : ^ aj8 (y) = ~^r9^(y) - ^i^g a u(y) (64') 
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o 



(65') 



For the same r-parameter family of infinitesimal diffeomorphisms of space- 



time manifold M onto itself {A(e 1 , . . . , e r ) 
we have 



A : M ->■ 
dy x 



| |e*| < 1,V1 < i < r}, 
6x K = ^-bx K (66) 

OX K OX K 

From this relation it is easy to see that Sx x (x; e 1 , . . . , e r ) and Sy x (y; e 1 , . . . , e r ), 
generally speaking, are functions of different forms. However, we will show 
that the continuity eqns. (65) and (65') are equivalent to each other. J£ and 
Jy are respectively the components in coordinate systems (x°, x 1 , x 2 , x 3 ) and 
(y°, y l ,y 2 , y 3 ) of the same vector field. 



Theorem 4 



d 



d 



dx 



(67) 



Proof. The lhs of eqn.(65) can be written as 

ux" uy 

v- 



dx 
dy x d 



dy 



dx 



TK 
•J 



dx K dy x 
d T / dv x 

dy 





dy 




dx 




dy 




dx 



d 



dx K 



dy 



dx 



d dy x 
dy x dx K 



V z \gW\J* 



dy 



dx 



+ 



cry dx a 
dx K dx a dy x 



dy } 



dx 
■dy x 



dx K x 





dy 




dx 



here an identity on the Jacobian (For its proof see appendix C) 



d 

dx K 



dy_ 

dx 



d f dy 



dy x \dx K 



dy_ 

dx 



(68) 



(69) 
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has been used. Hence the proof is reduced to proving 



In fact, it's easy to show that 5x K is a vector, 5u®(x) = u®(x) — u®(x) = 
5u e ^(x) — 5x p d p u^(x) is a (1, 1) -tensor, and 8g a p(x) = g a p(x) — Qapix) = 
Sg aj 3(x) — 5x p dpg aj 3(x) is a (0,2)-tensor. Hence the terms in the first brace 
of J£ eqn.(38) are vectors, and the first term in the second brace is a vector 
too. The rest of the terms in the second brace are not vectors individually. 
However, their sum is a vector. This can be proven straightforwardly, though 
tediously. (See Appendix B) ■ 

Conclusion 5 Now we have proven: 

(i) It is not the Noether's conservation laws written in different coordinate 
systems corresponding to the infinitesimal coordinate transformations with 
the same form that are equivalent to one another; rather, it is the Noether's 
conservation laws corresponding to the same family of infinitesimal diffeo- 
morphisms of spacetime onto itself that are equivalent to one another. 

(ii) All the Noether's conservative currents are vector fields (eqn.(70)), 
which should be the density, and current density of some scalar. Therefore, 

(Hi) All the Noether's conserved quatities are scalars. 

It is worth noting that the form of function J K eqn.(38) as a function 
of u(x), du(x), g(x), dg(x), d 2 g(x), Sx,Su(x),Sg(x),dSg(x) is independent of 
coordinate systems. However, for a given family of infinitesimal diffeomor- 
phisms of spacetime, the function forms of 5x(x, e 1 , . . . , e r ), 5u(x, e 1 , . . . , e r ), 
5g_(x, e 1 , . . . , e r ), d~5g(x, e 1 , . . . , e r ) and 6y(y, e 1 , . . . , e r ), Su(y, e 1 , . . . , e r ), 8g(y, e 1 , 
d5g(y, e 1 , . . . , e r ) are generally different, hence the function forms of J% and 
Jy as functions of field quantities, their derivatives, and coordinates, are dif- 
ferent. This does not ruin the covariance of our presentation which is assured 
by the coordinate independence of eqn.(38). 

We have proven these for the case of GR.. However it is easy to show 
they are true for the cases of SR and pre-relativity dynamics. 

According to Einstein, "[w]hat we call physics comprises that group of 
natural sciences which base their concepts on measurements; and whose con- 
cepts and propositions lend themselves to mathematical formulation." The 
aim of all natural sciences is to search for the objective laws of nature. 
In order to measure physical quantities or to describe physical processes in 
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terms of mathematics, one needs first to choose a reference coordinate sys- 
tem which depends on the observer's subjective will. Therefore, the ways of 
performing measurements and the formulations of physical laws, expressed 
in all reference coordinate systems should take the same form. This idea 
has guided Einstein from pre-relativity physics in Newtonian absolute space- 
time, to SR in Minkowski spacetime, and finally to GR in curved space- 
time. In particular, to get the proper motion equation one has to use an 
invariant Lagrangian C(g(x), dg(x), d 2 g(x),u(x), du(x)),such that the action 
A[g,u] = f x{n) d^x^J '-\g{x)\C{g{x), dg(x), d 2 g(x),u(x), du(x)) is an invari- 
ant under general transformation of coordinates. The least action principle 
would lead to a motion equation independent of coordinate systems. The gen- 
eral covariance of a physical theory means the operations of measurement, 
and the mathematical expressions of physical laws should be independent 
of reference coordinate systems. It does not mean that functions of coordi- 
nates, field quantities, and their derivatives, in different coordinate systems 
but with the same function form, must represent the same physical geomet- 
rical object, as is generally accepted. This is not always true as has been 
shown above. So, Noether's conserved current (38) 

J K [u(x),du(x),g(x),dg(x),d 2 g(x); 5x,5u(x), 5g(x),dSg(x)~\ =: J£ and 

J K [u(y),du(y),g(y),dg(y),d 2 g(y); 5y,Su(y),5g(y),dSg(y)] =: J* are com- 
ponents of the same vector field on spacetime; while the canonical energy- 
momentum (41) Tp(x) = Tp(u(x),du(x),g(x),dg(x),d 2 g(x)) and r K p {y) = 
Tp(u(y),du(y),g(y),dg(y),d g(y)) are not components of the same tensor 
field on spacetime. The long standing pseudotensor problems like (61), are 
only a result of misreading Noether's theorem; and the non-locality difficulty 
of gravitational energy-momentum does not really exist at all, had we read 
Noether's theorem properly. They are results of mistaking different physical 
geometrical objects as one and the same. 

It is important to distinguish a general law of nature and a concrete phys- 
ical object or process. The former is independent of coordinates, while the 
latter might looks different to different observers (from different coordinate 
systems). Yet we still can observe and describe a concrete instance in a way 
independent of coordinate systems. 

Now we are in a position to address the problems of conservation of 
energy-momentum and the gravitational energy-momentum in GR. These 
will be done in a later paper. 
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A Appendix 

The Lagragian density of classical fields L, is a function of the coordinates, 
field quatities, and their derivatives of up to the iV-th order. However, be- 
cause not all the variables are independent, such as d^d u u^(x) = d u d^u^(x) , 
9a/3{%) — 9pa{%), etc., there are infinitely many different function forms for 
L. This causes indefiniteness of derivatives, such as — 3 f/ - L, a 3 , If 

we drop all the redundent variables, then the Einstein summation convention 
can no longer be used, and the expressions will become awfully complicated, 
especially for a large N. In order to keep the formulaes neat, physists usually 
treat it in a different way. Here we will illustrate their method by using the 
lagrangian density for vacuum Einstein's equation, R (Ricci's scalar curva- 
ture). 

R is a function of 16 g a p(xys, 64 d fl g a p(xys, and 256 d IJ ,d ll g a p(xys. Be- 
cause g a p{x) = g/3a(x), and d^d^ga^x) = d u d^g a/ 3(x), there are only 150 in- 
dependent variables among them. We will choose 10 g a p(xys, 40 c^<7 a/ 3(x)'s, 
and 100 d lx d l/ g a ^(x) , s (a < f3,fx < v), for the independent variables. As a 
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function of 336 variables (As a function defined on a 336-demensional do- 
main), R can take infinite different forms, say, tp,ip,. ■ .When restricted to 
the 150-dimensional "sub-domain" D, all of them are the same function of 
150 variables. 

R\d = <p\d = ip\ D = ■ ■ ■ (Al) 

Substituting the 150 independent variables for all the variables in ip,ip, . . ., 
we get a unique function 

R{g a p{x),d^g al 3(x),d^d v g a p(x)), (a < (3 , fj, < v) (A2) 

Substituting \{g a p{x) + g^x)), ^(d li g a p(x) + d li gp a (x)), and ^d^g^x) + 
d^d v gp a {x) + d u d^g a p(x) + d u d^gp a (x)) for g a p{x), d^g a/3 (x), and d^d v g a ^(x) 
in R, respectively, we get a unique function of all 336 variables, denoted 
by R(g(x),dg(x),d 2 g(x)). This "standard" R(g(x),dg(x),d 2 g(x)) has the 
following property. 

dR 

dg a0 (x) 

dR 
dg aa (x) 



dR 



dR 



dgp a (x) 

dR 
dg aa (x) 



2dg a p{x) 



a < P 



(A3) 



dR 


dR 


dd^g a0 (x) 


dd^gp a (x) 


dR 


dR 


dd^g aa (x) 


ddf,g aa (x) 


dR 


dR 


dd fl d u g a0 (x) 


ddnd v gp a (x 


dR 


1 dR 


ddyd^g^x) 


4:dd^g a p{x) 



dR 



dR 



-,a < (3 



dR 



dR 



a < p, ii < v 
1 dR 



dd^gaaix) 

dR 
dd^ t d fM g a0 (x) 

dR 

dd tJi d IJi g aa (x) 



dd^gaaix) 

dR 
dd^gapix) 

dR 

ddpdngaaix) 
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2dd lx d 1/ g aa {x) 

1 6R 

2 ddpdpgapix) 



,a < (3 



(A4) 



(A5) 



(A6) 



When calculating the derivatives of R, we pretend that all its 336 variables 
are independent. Thus the indefiniteness problem no longer exists. 
From (Al), we have 

Slp\ d = 6ip\ D 

While 

dip dip dip 

r v-^ dip >r-^ d(f dip . . , . 

/x / dip 



a dd^gaaix) ^a</3 dd lx g a p{x) 

— 7 ; )Sd»gap(x) + V — ^ r Sd^gaaix) 

dd^gp^x) ^a,n dd^gaaix) 

dip dip 
^a<fs,^ dd^g a p(x) dd^d^gp a {x) 
dip dip 
+ ^ a ,ix<^dd^d u g aa (x) + dd u d^ 1 g aa ( y x)' >Sd ' ldu9aa(yX " > 

+ y ( *t + *e. + 

^ a <a^<v dd^g^x) ddvdpgapix) 

+ , , )5d,d u g aP {x)]\ D (A7) 



dd,j,d v gp a {x) dd u d^gp a (x) 

Because all the variations on the RHS of (A7) are independent, we get 

( dy? + dip = dj> + djj 
dg a p(x) dgp a (x) ° dg a ^(x) dgp a (x) ° 

( d( p + ft? )i = ( w + gV; )i 

dd^g^x) dd^g Pa (x) D dd^g^x) dd^g^x) D 

dd^g^x) ddvd^gapix) dd^g^x) dd v d^g Pa (x) D 

29 



^ dd p d u g aj3 (x) dd v d p g aj3 (x) dd^g^x) dd^g^x) ° 



This tells us, say, 



dg aj3 (x) 11 dg aj3 (x) dg a/3 (x) 

-5d^g a ^{x))\ D = ( — — ^ ^ a<3 (x))| D 



^dd p g aP (x) dd p g a/3 (x) 

#y^,0 a/ 3(x))|£> = ( a» ^ 7 ^ dpdvgapix)) 



^dd^d v g aP {x) ddnd v g a p(x 

dtp 



8dpd v g aP {x))\ D (A9) 



and even more, such as 



dd K g aP (x 
dip 



-){8g aP {x) - 5x p d p g aP {x))\ D 
■){5g a p(x) - 5x p d p g a p(x))\ D 



dd K g aP {x) 

where R(g(x), dg(x), d 2 g(x)) is the "standard" expression for R, and <p(g(x), dg(x), d 2 g(x)) 
is any expression from (Al). 

B Appendix 

Proposition 7 Let 

x ' [ dd K g a p(x) "dd K d,g aP (x) ^ ) dd K d p g^(x) ) 

dR 

(Sg a p(x) - 8x p d p g a p{x)) + — - - d^Sg^x) - 5x p d p g al3 (x)). (Bl) 



Then 

' l ~ i).r" 
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= ^ (B2) 



Proof. 



3 11 

-^d a g^dpg ar) - -d^g a/3 d v g pa - -d p g a( d v gf3 a - d a gp p d^g ar) ) (B3) 
f)R 1 

" ! = ^9^ + 9^9^) (B4) 



dd K d^g al3 2 

Note that g 9 g fl g is a (4, 0)-tensor, symmetrical for and for (a,/3). 

f ' = [g aK g^ + g^g<g^ + - <7<V V" 



-\g<*P g Wgt,v _ g< - g aR g^g m - g a ^g pK g m ] 
= :d lx g^B Ka ^ (B5) 



1 

OOKOfj.ga/3 Z 

+^(g aK d^ + g^d^) - g a ^ 



+d,ga-\9 aK 9^9 m - \<f*9 K <r + 9^9^9 m ] 
= : d^A"*™" (B6) 

Note that A Ka ^ and B Ka ^v are (6, 0)-tensors.Let 



= -\g a ^x)g^{x)g^{x) + ±<r(x)gK(x)gr(x) 
-\sr*(x)tf*{x)gr{x) - \g m {x)s^{x)sr{x) 

+ \g a Kx)9^{x)g^{x) + g aK (x)g^(x)g^(x) - \g af) {x)g^{x)g^{x) (B7) 
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Then 



0*) = t^ T (,) + A(gl) (B8) 



d^g^x) 



i OR d dx^ dR 

~ K '^ y) dd^g a ,^y) + ^W^O^g^iy/ 9 ^ 

dx K dR dx a dx 13 d dx a dx 13 - 

— OR — Ox K 

dx K r „i dx K 

= ^r^K + ^^ rest 
ay dy K 

d . dx^ dR - 

dx» dy» ddn'd^g^piy) 

We are going to show that rest vanishes. Its first term is 

= [g a ' e (y)g"' K ' (y)g" w (y) - \g a ' e (y)g" w (y)g K '"' (y) 
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where 



-g a ' K/ (y)g^'(y)g^'(y) + \g a ' p '{y)g^\y)g^\y)] * 

d 2 x^ dy^ 



gt"v>(y) s ga'p'(y) 



= g a ''\y)g"' K '{y)^{^)~8gMy) (A) 

1 f) 2 T^ Bill 3 ' - 



-g^\y)g^\y)^{^gMv) (A) 
+\^(v)g""(v)^(^)hMv) ((C) (bio) 



The second term is 



= (y)/ v (y)/"' (y) - ^ (y)/"' (y) <?" v (y) 



-^(v^W^i/) + ^ Q ' /3 '(y)/ /1 '(y)/' ? '(y)] x 



dy^'dy v ' dx v 



gev"(y) s 9a'^'(y) 



g"' K \y)g^\y)Q^%fhMy) (d) 
-^^'(y)^'(y)^^^W(y) (B) 



-y a ' K/ (y)/ v (y)^(^)^ r (y) (A) 
+^ a,/3 '(y)y^'(y)^(^)W(2/) (c) (bh) 
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The third term is 



d , dx^ dR — 



%){9 a ' K 'g^'-9 a '"'g^¥9My) 



= ^w )ga ' K,/ "~ 69M ^ (A) 

( c ) ( Bi2 ) 



The forth term is 



dx^ dyv' 
OR d 2 x a dy a ' - 



9 a " K '9^'5gMy) (B) 



dy^'dy a " dx c 

+ ww^^' TS9My) (B) (B13) 



The last term is 

dR d 2 x? dyP' 



Sga/p/(y) 



dd K >d li ig a >p«(y)dyi i 'dyP a dx? 

(gartgfrS -g«'P"g«'»')8g a ,p,(y) 



d 2 x 13 dy 13 ' , lR , _ /( g„ k'h' 



dy^'dyP" dx$ 

''-.'/""'// ; " %,•,•!//) (D) 



ch/' 3 " cte' 3 ' 

+ WW^ S'"S9a'Ay) (B) (B14) 

All the terms marked (A) cancel each other, all the terms marked (B) cancel 
each other, etc. Therfore we get 

4 K = |k' (B15) 

That is (B2). ■ 
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C Appendix 



Proposition 8 



d 

dx K 



dy_ 

dx 



d ( dy 



dy x \dx K 



dy_ 

dx 



Proof. 



d 
dx K 



dy_ 

dx 



d dy } 



d 



dx K dx i 
d 2 y x dx' 



Pi( dy x 



dx° 

dy 



dx K dx a dy } 
d dy 



dx 



dy x ^ dx K 



dy_ 

dx 



dy_ 

dx 
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